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Abstract Gamma-rays and neutrons are the only sources of information on energetic ions 
present during solar flares and on properties of these ions when they interact in the solar 
atmosphere. The production of y-rays and neutrons results from convolution of the nuclear 
cross-sections with the ion distribution functions in the atmosphere. The observed y-ray and 
neutron fluxes thus provide useful diagnostics for the properties of energetic ions, yielding 
strong constraints on acceleration mechanisms as well as properties of the interaction sites. 
The problem of ion transport between the accelerating and interaction sites must also be 
addressed to infer as much information as possible on the properties of the primary ion 
accelerator. In the last couple of decades, both theoretical and observational developments 
have led to substantial progress in understanding the origin of solar y-rays and neutrons. 
This chapter reviews recent developments in the study of solar y-rays and of solar neutrons 
at the time of the RHESSI era. The unprecedented quality of the RHESSI data reveals y- 
ray line shapes for the first time and provides y-ray images. Our previous understanding 
of the properties of energetic ions based on measurements from the former solar cycles is 
also summarized. The new results - obtained owing both to the gain in spectral resolution 
(both with RHESSI and with the non solar-dedicated INTEGRALISPl instrument) and to 
the pioneering imaging technique in the y-ray domain - are presented in the context of this 
previous knowledge. Still open questions are emphasized in the last section of the chapter 
and future perspectives on this field are briefly discussed. 
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1 Introduction 

The first evidence of acceleration of charged particles (electrons and ions) to relativistic en- 
ergies in association with solar flares was found in 1942 with the discovery of large increases 
in the count rates of several ground level cosmic-ray intensity monitors. These observations 
told us, more than fi fty years ago, th a t sola r flares are capable of accelerating protons to 
GeV energies. While iBiermann et alJ dl95lh had predicted that relativistic protons accel- 
erated during a solar flare could produc e a flux of high -energy neutrons observable at the 
Earth (as well as some y-ray emission). iMorrisonl |l958) published the first prediction of y- 
ray fluxes expected from several celestial objects including solar flares. These first estimates 
were at the base of for mer detailed theoretical studies of high-energy neutral emissions 
which were initiated by Lingenfelter & Ramatvl dl967ah . They provided a detailed treat- 
ment of the production of nuclear reactions in the solar atmosphere and calculated the yield 
of neutrons, positrons and y-ray lines using normalized rigidity spectra for the accelerated 
ions. The first observational evidence for y-rays cam e from the observa tions of the flare 
SOL1972-08-04 with the OSO-%\ y-ray spectrometer dChupp et al.lll973l . see Figure [TTTI 
which confirmed the prediction by Morrison that nuclear reactions in the solar atmosphere 
could produce y-ray lines. Pioneering observations of solar y-ray lines were also obtained 
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Fig. 1.1 Time-integrated y-ray spectrum observed aboard OSO-7 in the solar and background quadrants for 
the discovery event, SOL 1972-08-04. The integration time is from 06:24 to 06:33 UT. The solar quadrant 
shows significant enhancement up to 7 MeV as well as line emissions at 0.5 MeV, 2.2 MeV, 4.4 MeV and 
6.1 MeV. The lines at 1.17 and 1.33 MeV are from the on-board calibration source (from Suri et al. 1973). 



for later flares of August 1972 by spectrom eters on Prognoz-6 llTalon etal.ll975h and for the 
flares SOL1 978-07-11T10:54 b y HEAO 1 faudson et alJll980h and SOL1979-11-09T03:20 
by HEAO 3 jPrince et alj|l982h . Since then solar y-ray astronomy has provided a new win- 
dow for the study of ion acceleration in solar flares. 

Energetic ions interacting with the solar atmosphere produce a wealth of y-ray emis- 
sions. A complete y-ray line spectrum is produced through interactions of ions in the range 
~1-100 MeV/nucleon and consists of several nucl ear de- e xcitati o n lines, neutron capture 
and positron annihilation lines. See, for example, Chupp (1984), Ramatv ( 1986), Chupp 
jl996h.lTrottet & yilmerl l ll997l) . IShare & Murphvl ( l200fj|)lvilmer & MacKinnon! j2003h . and 
IShare & Murphvl d2006h . The temporal and spectral characteristics of all these radiations 
provide strong constraints on acceleration timescales, ion energy spectra and numbers as 
well as energetic ion abundances. If more energetic ions are present, over a few hundred 
MeV/nucleon, nuclear interactions with the ambient medium produce secondary pions whose 
decay products lead to a broad-band continuum at photon e nergies above 10 MeV (with 
a bro ad peak around 70 MeV from neutral pion radiation) (Stecker 1970; Murp hy et all 
Il987l) and also secondary energetic neutrons which, if energetic enough, may escape from 
the Sun and be directly detected in interplanetary space (> 10 MeV neutrons at 1 AU) or 
at ground level (> 200 MeV neutrons) (see, e.g., IChupp & DunphvlbOOdJchupp & Rvad 
200 9). The importance of obs erving neutrons was pointed out bv lBiermann et alJ 1 1951 ) and 
by iLingenfelter et al.1 Ill965al lbh. Lower-energy neutrons (kinetic energies below lOOMeV) 
can only be observed in space because they are strongly attenuated in the Earth's atmo- 
sphere and cannot reach the ground. Neutrons with kinetic energies higher than 100 MeV 
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can be observed on the ground and hence with simultaneous observations in space and on 
the ground, it is possible to obtain the energy spectrum of solar neutrons and of accelerated 
particles, in a wide energy range. 

The first detection of solar high-energy neutrons was provided by the observations of 
the 5MM/GRS0 spectrometer for the flare SOL1980-06-21T02:00 (X2.6) which showed 
not only the emission from ultra-relativistic electrons at the ti me of the flare but also 
a later signature from high-energy neutrons dChupp et al.l 1 1982). This same experiment 
provided the first observations of pion-d ecay radiation in s olar flares (high energy radi- 
ation produced by ultra-relativistic ions) dForrest et al] 1 19851) for the event SOL 1982-06- 
03T13:26 (X8.0). For this event, solar neutrons were also detected for the first time at 
ground-level by the I GY-type neutron m onitor installed at Jungfraujoch, Switzerland, as 
well as by SMM/GRS (Chupp et al. 1987). With these observations, all t he component s pre- 
dicted by Lingenfelter & Ramatv ( 1967a) had been detected (see, e.g.. lChuppll996S . for a 
detailed review of early y-ray observations of solar flares). Since these first observations, 
solar energetic neutrons have been detected with ground-based detectors on a handful of oc- 
casions. 10-120 MeV neutrons were also detected in space with COMPTEL on the Comp- 
ton Gamma-Ray Observatory ( CGRO) mission, whose mode of oper ation enabled direct 
determination of their energies llRvan et al.ll 19931 : iKocharov et alJ fl998). We consider these 
various neutron observations in more detail below (Section[4]i. 

Before the launch of RHESSI and INTEGRAL in 2002, high energy observations of so- 
lar flares had been obtained for two solar cycles respectively with SMM/GRS as well as 
on Hinotori for cycle 21 and by the PHEBUS and SIGMA experiments aboard Granat, the 
wide-band spectrometer on Yohkoh, GAMMA-1 and CGRO for cycle 22 (see, e.g., Chupp, 
1984; Chupp, 1996; Ramaty, 1986; Trottet & Vilmer, 1997; Vilmer & MacKinnon, 2003, 
for reviews). Before the advent of RHESSI observations, the solar electromagnetic radiation 
above 100 keV (including the y-ray line region) was one of the last spectral domains with no 
spatially resolved observations, so that no information was available on the location of the 
interaction sites of the energetic ions (nor of electrons above 100 keV) on the Sun. y-ray line 
emission from ions had been observed in many flares but the quantitative constraints from 
y-ray line spectroscopy, which provides information on the bulk of the flare fast ions, whose 
presence had been deduced from observations with limited spectral resolution. Higher en- 
ergy radiations and neutrons have been observed in fewer events, but they provide crucial 
information on the extremes of particle acceleration. 

Since the laun ch of the solar-dedica ted RHESSI mission fcin et al.ll2002h and of IN- 
TEGRAL in 2002 dVedrenne et alj [2003), both equipped with high-resolution germanium 
detectors, several y-ray line flares produced strong enough fluxes so that a detailed y-ray 
line shape analysis could be achieved for a few events. Furthermore, since February 2002, 
RHESSI h as observed a tota l of 29 events above 300 keV with 1 8 clearly showing y-ray line 
emission dShih et alj|2 009b). Pioneering results were obtained by RHESSI on a few of these 
events on the localization of the interaction sites of energetic ions. We shall review below 
the results provided by RHESSI and also by INTEGRAL on these issues. 
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Fig. 2.1 OSSE spectrum of the flare SOL1991-06- 04T03:37 (X12.0) summ arizing the different components 
and the physics from y-ray line spectroscopy (from Share 

2 y-ray line observations of energetic ions 

2. 1 General features of a y-ray spectrum 

Many studies performed before the launch of RHESSI have been devoted to the study of 
the complete y-ray line spectrum which is produced in the solar atmosphere through the in- 
teractions of ions in the ~1-100 MeV/nucleon range. In contrast to the monotonic electron 
bremsstrahlung continuum, the nuclear y-ray spectrum displays a rich variety of features: 
broad and narrow lines and various sorts of continuum arising from several mechanisms. The 
bombardment of ambient nuclei of the solar atmosphere by accelerated protons and a parti- 
cles produces a whole spectrum of de-excitation y-ray lines (prompt narrow lines), almost all 
of them with rest energies lying between a few hundred keV and 8 MeV. Inverse reactions 
in which accelerated C, O and heavier nuclei collide on ambient H and He are the origin 
of prompt broad lines centered at the same line energy. Since the emitting nuclei retain 
significant momentum after the exciting collisions, these lines are substantially Doppler- 
broadened and/or shifted, to the extent that they are not individually resolvable. These broad 
lines merge with many other weaker lines from nuclei he avier than oxyge n to form an unre- 
solved continuum (see Section l2.1U ( Murp hy" et al.l20 091. Figure are & Murphy 

feOOCh . based on the HEAO-1 /OSS$ observations of the large flare SOL1991-06-04T03:37 
(X12.0), shows the major features of the y-ray line spectrum: strong narrow de-excitation 
lines are found at 6.129 MeV from 16 0, 4.438 MeV from 12 C, 1.779 MeV from 28 Si, 
1.634 MeV from 20 Ne, 1.369 MeV from 24 Mg and 0.847 MeV from 56 Fe. The broad lines 
merge into a quasi continuum dominating over the bremsstrahlung emission from electrons 
in the ~l-8 MeV range. The single strongest line is the extremely narrow 2.223 MeV deu- 
terium formation line. Also present are photons of flare origin at 0.5 1 1 MeV from positron 
annihilation and the broad feature at 0.4-0.5 MeV formed from the two lines, at 0.429 and 
0.47 MeV, that accompany formation of 7 Li and 7 Be in a- 4 He fusion reactions. 
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2.2 Magnetic loop transport and interaction model 

2.2.1 Gamma-ray line production including transport 

Information derived from y-ray line emissions on the characteristics of the accelerated ions 
requires a comprehensive understanding of the transport of these energetic ions between the 
acceleration site, presumably located in the low corona, and the interaction and emitting 
sites in the denser layers of the solar atmosphere. This requires a model for ion acceler- 
ation, transport in flare loops and interaction. Several models have been developed since 
the 1980s on this to pic, a couple of which we describe. The recent model described by 
Mur phy et alj d2007l) addesses particle transport and interaction and includes energy losses 
due to Coulomb collisions, removal by nuclear reactions, magnetic mirroring in the converg- 
ing magnetic flux and MHD pitch-angle scattering in the corona, tre ated via the quasi linear 
formalism. The loop formalism is the one previously developed bv lHua et al.l (fl989). The 
loop consists of a semicircular coronal portion of length L and two straight portions extend- 
ing vertically from the transition region through the chromosphere into the photosphere. 
Below the transiti on region, the magnetic field strength is assumed proportional to a power 
8 of the pressure dZweibel & Haberill983l) taken here as 5 ~ 0.2. Such a converging mag- 
netic field results in mirroring of accelerated particles. Pitch-angle scattering is characterized 
by A , the mean free path required for an arbitrary initial angular distribution to relax to an 
isotropic distri bution. The dep endence of A on particle energy is expected to be weak (see 
discussion bv lHua et al.ll 19891) and is assumed to be independent of particle energy. In the 
model, the level of pitch-angle scattering of the energetic particles is simply characterized 
by A, the ratio of A to the loop half-length L c = L/2. Magnetic convergence and pitch-angle 
scattering determine the angular distribution of the accelerated particles when they interact 
with the ambient medium, which is crucial to several aspects of the observable emissions. 
Several height-density profiles (n(h)) for the solar atmosphere are also assumed. 

2.2.2 Time history of the nuclear interaction rate 

iMurphv et all J2007I) studied the effects of transport on temporal behavior by assuming in- 
stantaneous release of all particles at t = 0. Actual observed time profiles represent the con- 
volution of this behavior with the time profile of the particle accelerator. In the absence of 
magnetic convergence, ions do not mirror and the nuclear interaction time history depends 
only on the energy loss rate in the lower chromosphere and upper photosphere where the 
density is greatest and most of the interactions occur. 

In the presence of magnetic convergence, particles behave differently depending on 
pitch-angle a. Increase of the magnetic field and of the density towards the deeper at- 
mosphere leads to partial particle trapping and to the formation of a loss cone of angular 
half-width ((Xq). Ions with a < Oq precipitate in a loop transit time, behaving as described 
above. Ions outside the loss cone mirror and lose energy much more slowly as they traverse 
the low-density corona. These interactions thus occur on longer time scales, with the time 
scale increasing with increasing convergence. Pitch-angle scattering causes the loss cone 
to be continuously re populated, and in fact is essential if obs erved, impulsive time profiles 
are to be understood jZweibel & Haberll 19831 : iHua et al.ll 19891) . As a result, with increasing 
pitch-angle scattering rate (i.e., decreasing A) the nuclear interaction rate increases at early 
times and correspondingly decreases at later times. However, the time history is no longer 
affected by increasing pitch-angle scattering when the rate of loss-cone replenishment ex- 
ceeds the rate of nuclear reactions in the loss cone. Increasing the loop length increases the 
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time scale of the interaction rate since mirroring particles spend more time at lower coronal 
densi ties where nuclear re actions are less likely. All these effects are illustrated in Figure l2~2l 
from lMurphv et alj [2007 ) . On short timescales (<1 s) the effects of velocity dispersion may 
be seen, as higher energy particles transit the loop and encounter the deeper regions more 
rapidly. On longer timescales, trapping and precipitation, via pitch-angle scattering, domi- 
nate the temporal behavior. 

The results of these detailed calculations are summarized and put to work in the "trap- 
plus-precipitation" model, in which the flaring atmosphere is divided into a coronal, mag- 
netic trap characterized by a single, uniform density no, and a high density, chromosphere- 
photosphere region. The escape of particles from the trap to the deep atmosphere is sum- 
marized in an escape rate. The saturation of escape rate with increasing scattering rate, 
de scribed above, may be rec ognized as the "strong pitch-angle scattering" limit identified 
bv lKennel & PetschelJ dl966h . in which case the coronal pitch-angle distributi on is always 
more or less isotropic and the escape time from the trap is given by (see, e.g.. lHulot et alj 
1 19891) 

/y2 y 

JW = y^ (2.1) 

where (Xq and L are as above and V is the particle velocity. For scattering rates below 
the strong limit, the resulting long coronal particle lifetimes become incon sistent with the 
short duration usually observed in solar y-ray bursts dZweibel & Ha ber 1983). The trap-plus- 
precipitation model has been applied to observations, with results described below (Section 

i27za . 

2.2.3 Angular distribution of the interacting ions 

In the presence of magnetic convergence, particles encounter the greatest densities at their 
mirror points when they also have the greatest pitch angles. Observations (e.g., of line 
shapes) reveal the density-weighted, source-in tegrated pitch-angle distribu tions and these 
may consequently take a "fan-beam" form (e.g jHua & Lingenfelterill987bT) . As pitch-angle 
scattering approaches the strong limit (A w 20), however, more and more particles precip- 
itate in the loss cone, at small pitch angles, and the density-weighted distribution becomes 
more strongly beamed downward. So the overall pitch-angle distribution reflects a combi- 
nation of the initial pitch-angle distribution, the degree of magnetic convergence in the loop, 
and the pitch-angle scattering rate (see Murphy et al. 2007 for details). 

2.2.4 Depth distribution of the interaction site 

In the absence of magnetic convergence, there is no mirroring and the depth distribution 
directly reflects the grammage required for the accelerated ions to interact as they move 
downward through the solar atmosphere. In the presence of magnetic convergence, mirror- 
ing results in interactions occurring at higher elevations (and therefore lower densities) as 
more and more particles are prevented from penetrating the lower atmosphere. Even with 
minimal convergence and in case of no pitch-angle scattering, a significant fraction of the 
interactions occurs at low densities (almost 20% of the interactions occur at densities less 
than 10 10 cm -3 ). As the scattering is increased, more particles are able to precipitate and 
the bulk of the interactions move deeper, to higher densities. Because higher-energy ions 
tend to interact farther along their paths, harder particle spectra result in interactions oc- 
curring at higher densities since more higher-energy ions are producing the interactions. 
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Fig. 2.2 Nuclear interaction rate as a function of the different parameters: (a) magnetic convergence, (b) 
pitch-angle scattering, (c) loop length, (d) spectral index, (e) atmopsheric model, (f) combined effect of the 
pitch-angle scattering and of the atmospheric model. The yields are normalized to one proton of energy 
greater than 30 MeV (from M urphy et all2007n . 



Similarly, interactions whose cross-sections have higher threshold energies tend to occur 
preferentially at higher densities. This dependence of the depth distribution on the inter- 
action cross-section also explains why the depth distribution is affected by the accelerated 
a/proton ratio. Alpha-particle interactions generally have lower threshold energies and, as 
just discussed, such interactions occur at lower densities. When the a/proton ratio is high, 
a larger fraction of the line yield is due to such reactions and the depth distribution shifts to 
lower densities (see lMurphv et alJl2007l . for details). 

2.2.5 Effective energies of the accelerated ions producing y-ray lines 



The main factor determining the effective energy range for a nuclear reaction is the cross- 
section, but the shape of the accelerated-ion spectrum can be important. For thin-target in- 
teractions, only ions with initial energies where the cross-section is significant contribute. 
For thick target interactions, however, even ions with higher initial energies also contribute, 
since they can lose energ y and then interact at energies where the cross-section is significant. 

iMurphv et all fc007t) determined the effective ion energies to produce a y-ray line for 
a particular ion spectrum by weighting yields from monoenergetic ions with the spectrum. 
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Fig. 2.3 Yield of the 6.129 MeV 16 line weighted by accelerated ion power-law spectra with indices s = 
3 and 5. The contr ibutions to the line from accelerated protons and a particles are separately shown (from 
iMurphv et all2007l) . 

The ambient abundances are assumed to be coronal teamatvetaT1ll996h . As an example, 
Figure [231 from IMurphv et all J2007I) shows the differential yield of the 1.634 MeV 20 Ne 
line as a function of accelerated ion energy for power-law spectra with indices of 3 and 5 and 
accelerated a/proton = 0.5. These energies are the energies of the ions leaving the accelera- 
tion region. Separate contributions to the line from accelerated protons (dashed curves) and 
a particles (dotted curves) are plotted. For soft spectra, the a-particle contributions domi- 
nate due to their lower threshold energies, and the most effective ion energies are around a 
few MeV/nucleon. For hard spectra, the proton interactions become important due to their 
higher threshold energies and because the effective energies are higher. An effective ion en- 
ergy range for producing the line is defined as follows: the ion energy where the yield is 
maximum is determined and the effective range is the one for which the yield has fallen to 
50% on each side of the maximum. We note that if the effective energy distribution is very 
broad, the arbitrary value of 50% could be misleading. 

The fluence ratio of the 6.129 MeV 16 and 1.634 MeV 20 Ne lines has been frequently 
used as a measure of the accelerated-ion spectral index. Figure 12.41 shows the effective 
accelerated-ion energy ranges for producing the 1.634 MeV 20 Ne line (white boxes) and the 
6.129 MeV 16 line (grey boxes) for several power-law spectral indices and for accelerated 
a/proton = 0.1 (panel a) and 0.5 (panel b). The horizontal line within each box corresponds 
to the peak of the distribution. The Figure shows that the effective energy range for produc- 
ing the 16 line is usually shifted to energies higher than for the 20 Ne line. The extent of the 
separation of the effective energy ranges for the two lines determines the sensitivity of the 
ratio to the spectral index. Because the separation for these two lines is not large, the ratio 
is not very sensitive. For very hard spectra, the upper range of the effective ion energies can 
be ~100 MeV/nucleon, but for most (softer) spectra it is less than 10 MeV/nucleon. The 
lower range is typically a few MeV/nucleon but for very soft spectra (particularly when the 
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Fig. 2.4 Effective accelerated ion energies for production of the 1.634 MeV 20 Ne line (white boxes) and 
6.129 MeV ie O line (grey boxes) as a function of spectral index. The effective energy ranges are the 50% 
yield range as defined in the text. Panel (a) is for accelerated a/proton = 0.1 and panel (b) is for 0.5 (from 
Murphy et al. 2007). 



a/ p ratio has the commonly assumed value of 0.5) it can be less than 2 MeV/nucleon. As 
the spectrum hardens, the higher-energy proton reactions begin to contribute and the effec- 
tive ion energy range shifts to higher energies and becomes broader as interactions of both 
a particles and protons contribute. The interpretation of this particular energy-distribution 
diagnostic is discussed further below (Section [23}, where the roles of ambient chemical 
abundances and the relative contributions of protons and a -particles are discussed. 

The fluence ratio of the 4.438 MeV 12 C de-excitation line and the 2.223 MeV line has 
also been used as a measure of the accelerated-ion spectral index. Figure l2~31 shows the effec- 
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tive accelerated ion energy ranges for the 4.44 MeV line (white boxes) and the 2.223 MeV 
line (grey boxes) for several power-law spectral indices and for accelerated a/proton = 0. 1 
(panel a) and 0.5 (panel b). For hard spectra (spectral index s <4), the effective ion energies 
for the neutron capture line extend to very high energies, up to ~ 100 MeV/nucleon, and are 
generally much higher than those producing the 12 C line. However, for soft spectra (.? > 4), 
the effective ion energies for the neutron capture line can be very low (<7 MeV/nucleon) 
due to neutron production by the exothermic (a,w) reactions on heavy elements. For such 
soft spectra the effective ion energies for the neutron capture line are even less than those 
for the 12 C line. As the flare location moves from disk center to limb, the effective energies 
shift to slightly lower energies since the neutrons from higher energy reactions are gener- 
ally produced deeper and are subsequently more attenuated by Compton scattering. When a 
spectral index is derived using the ratio of these two lines, the relevant ion energies cover a 
broad range of energies. For hard spectra, the relevant ion energy range extends from a few 
MeV/nucleon up to and greater than 100 MeV/nucleon. For soft spectra (s >4), the relevant 
ion energy range is much narrower, from around 1 to a few MeV/nucleon. For very soft 
spectra, the lack of separation of the effective ion energies producing the two lines reduces 
the sensitivity of the ratio to the spectral index. 



2.2.6 Loop parameters and time development of J-radiation 

iMurphv et ai] d2007l) discuss an integrated approach to the interpretation of y-ray measure- 
ments, distinguishing in particular between the deduction of quantities describing the accel- 
erated particles ("acceleration parameters") and quantities describing the atmospheric en- 
vironment in which the particles evolve and radiate ("physical parameters" - loop lengths, 
target abundances, parameters describing the state of plasma turbulence in the corona), y- 
ray measurements in turn yield various quantities, each of which contains information on 
some combination of acceleration and/or physical parameters. De-excitation line fluences, 
for example, are determined by acceleration parameters and otherwise only by source region 
chemical abundances. Line widths and centroid energies, on the other hand, are determined 
by the energy distribution and composition of energetic ions, but also by the amount of 
magnetic convergence in the loop and the degree of pitch-angle scattering, i.e., in princi- 
ple, the energy distribution over wavenumbers of coronal MHD turbulence, but typically 
parametrized in terms of the ion mean free path A (normalized to the loop half-length). 
Below we concentrate on deduction of loop parameters from timing measurements. 

Peak time delays, i.e., differences between the times of peak flux measured at different 
photon energie s, have been reported in a few fla res between hard X-rays and y-ray line 
emissions (e.g. jHulot et al.ll 1 9891 : [^"shimorill 1989h . Their interpretation was focused on the 
deduction of loop parameters. The relative timing of prompt y-ray line emissions and of hard 
X-rays was studied in this context by Hulot and collaborators. These models were again used 
to interpret the delays of the hard X-ray an d y-ray line emissions observed for two flares by 
ISEE-3 and SMM/GRS (Hulot et al.ll 19921) and a lso applied to one of the y-ray line flares 



observed with RHESSI (Dauphin & Vilmerl l2007l) . 

For the study of temporal behavior, the basic characteristics of the hard X-ray and y-ray 
trap- plus-precip i tation models are summarized above, and described in lVilmer et al.l ( 119821) 
and lHulotet^i] | |l989l) . The principle of the model is the following: electrons and ions are 
supposed to be continuously injected over a finite injection time in a coronal loop where they 
experience trapping and pitch-angle scattering as above. While particles trapped in the loops 
produce X-ray and y-ray thin-target emissions, electrons and protons escaping from the trap 
propagate downwards to the denser chromospheric layers where they produce hard X-ray 
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Fig. 2.5 Effective accelerated ion energies for production of the 4.438 MeV C line (white boxes) and 
2.223 MeV neutron capture line (grey boxes) as a function of spectral index. The effective energy ranges are 
the 50% yield range as defined in the text. Panel a is for accelerated a/proton = 0. 1 and panel b is for 0.5. For 
the neutron capture line, X = 300, 8 = 0.2, L = 1 1,500 km, and ob8 = 75° (from Murphy et al. 2007). 



and y-ray line emissions in a thick-target approximation. In these works, the precipitation 
rate is supposed to be in the "strong diffusion limit" given in Equation I l2.lt . The free, 
"physical" parameters of the model describing the loop thus reduce to the precipitation rate 
(Section r2.lt of the trapped particles into the chromosphere and the density of the ambient 
medium. Additional parameters describe the production of energetic ions: the spectral index 
and the time evolution of the injected particles which in this model are not instantaneously, 
but continuously released. Energetic electrons, protons and alphas are assumed injected in 



Gamma Rays and Neutrons 



13 




Time (sec) 



Fig. 2.6 Time profile of the production of the line at 4.4 MeV in the case of the trap-plus-precipitation model 
developed by Hulot et al. (1989) in the case of a short injection time t mcu and for two values of the precipitation 
rate ag/L = 1(T 5 km" 1 and KT 6 km" 1 (from Hulot et al. 1989). 



the corona at a rate given by 

q(E,t) = q f(E)f{t), (2.2) 

where qo is the amplitude of the injection and f(t) is given by a simple parametric form 
with rise-and-fall behavior: 

/(f) = (f - f )(2f mm - - - to)) for t < t < 2t max + t ; 
/(f) = elsewhere, 

where t max and to are, respectively, the peak time and the starting time of the injection. Here 
f{E) represents the energy spectrum of injected electrons and ions; f(E) is taken as a power 
law in energy with a spectral parameter 8 defined by 

f(E)=E~ s . (2.4) 

These models compute electron bremsstrahlung radiation and the strong y-ray lines at 
6.129 MeV from 16 and 4.438 MeV from the first excited state of 12 C, populated directly, 
and by the spallation reactions on 16 0, and 1.779 MeV from 28 Si, 1.634 MeV from 20 Ne, 
1.369 MeV from 24 Mg (with the photospheric target abundance for y-ray lines). The free 
parameters are finally the parameter a, 2 jh linked to the precipitation rate, the density, the 
starting and the peak time t max of the injection, the amplitude of the injection and the spec- 
tral index of the energy spectrum of injected particles. Figure 12.61 shows the effect of the 
precipitation rate on the delay with respect to the injection time and on the time decay of 
the production of the prompt y-ray line. Decreasing the value of the precipitation rate has 
a similar effect on the decay time as decreasing the pitch-angle scattering in the previous 
model. 

lHulotetaD fl989). by using a trap-plus-precipitation model, were the first to calculate 
the relative time delay of y-ray flux with respect to the hard X-ray flux. They showed that 
such delays strongly depend on both the trapping time of the accelerated particles and the 
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Fig. 2.7 Time evolution of the bremsstrahlung flux at 150 keV and 1 MeV and of the prompt nuclear de- 
excitation lines observed with a time integ ration of 20 s with RHE SSI for SOL20 02-07-23Tq0:35 (X 4.8). 
The starting time of the figure is 00:26 UT jDauphi n & Vilmer 2007), adapted from lShare et"aT] <2003all 



density of the loop. In particular, they showed that the time delays of the hard X-ray and y- 
ray line emissions observed by ISEE-^ and SMM/GRS during the SOL1980-06-07T03:22 
and SOL1981-04-27T09:45 events could be explained by this model and that parameters 
of the particles and of the ambient model could be obtained. However, these studies were 
performed using the time profile of the 4-8 MeV range as a proxy for the time profiles of the 
nuclear lines, and no images at X-rays or y-rays were used as additional constraints on the 
model. 

A similar analysis was performed bv lDauphin& VilmeJ d2007h on the first y-ray line 
event observed by RHESSI, SOL20 02-07-23T00:35 (X4.8). Indeed, the analysis of the hard 
X-ray and y-ray flux time evolution dlin et al .120031 : IShareet al.l2003ah showed that there is 
a slight delay between hard X-ray emission at 150keV and y-ray line emission. Figure |2~T 
shows the temporal evolution of the HXR emissions at 150 keV and 1 MeV and of the 
y-ray line ti me profile. These ti me profiles are obtained from the spectroscopic analysis 
presented in IShare et alJ d2003ah . The spectral analysis is performed for successive time 
intervals of 20 s, and results in the fitting to observed count spectra of a model photon 
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Fig. 2.8 Time evolution of the observed and com puted y-ray flux derived for different parameters of the trap 
plus precipitation model (from Dauphin & Vilmei; l2007l) . 



spectrum including a double power law for the bremsstra hlung continuum, a nuclear 7- 
ray line func tion made of 15 nar row and broad Gaussians ( Smi th et al.ll2003r) . a neutron- 
capture line (Murphv et al ] l2003h . an a— 4 He fusion line complex between ~400 and 500 
keV dShareetal .112 003a), and the solar annihilatio n line a t 511 keV and its positronium 
continuum. Figure 12.71 - adapted from IShare et al.l d2003al) - shows the results of this fit 
as a function of time: the bremsstrahlung flux at 150 keV and 1 MeV and the nuclear de- 
excitation line flux. This figure shows that the hard X-ray and y-ray time evolutions are 
roughly similar, indicating a common origin of the accelerated electrons and ions, but that 
for the first main peak (from s to 400 s) a time delay of around 12 s is observed between 
the time profile at 150 keV and the y-ray line time profile. There is, however, no significant 
delay between the time profile at 150 keV and at 1 MeV. 

iDauphin & VilmeJ J2007h investigated whether this delay could be reproduced by the 
trap-plus-precipitation model, given the additional constraint provided by the X-ray and 7- 
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ray images. They found that the time profiles of both HXR and prompt y-ray line fluxes 
could be reproduced (given that electrons and ions are injected and partially trapped in 
different coronal loop systems with slightly different characteristics discussed below) but 
they also found that the detailed X-ray and y-ray time profiles could be well reproduced 
if the ratio of energetic electrons with respect to ions slightly varies from peak to peak 
during the flare. Such a variability of electron-to-ion production from peak to peak had 
been observed previously in flares based on spectral analysis (e.g., IChupp et al.lll993h and 
iDauphin & Vil mer ( 2007) showed here that a similar result can be found from time-delay 
analysis. They finally found, given this variability in the electron-to-ion ratio, that a good 
reproduction of the evolution of both X-ray and y-ray line time profiles (see Figure l2~8t was 
obtained if compared to electrons - ions were injected in a system of coronal loops of lower 
density and with a slightly larger mirror ratio leading to a slightly more efficient trapping. 
If we consider a similar value for the loss cone for the two coronal loops, we find that 
ions should propagate in slightly longer loop lengths, which is consistent with the imaging 
observations (see Section[3]l. 

The emphasis abo ve has been on inter-c omparison of X-ray and y-ray time profiles. As 
discussed in detail bv lMurphvetal.N2007h . analysis of y-ray measurements on their own 
can also provide constraints on the loop parameters, e.g., from the flare y-ray line shapes 
and the fluence and time evolution of the 2.223 MeV neutron-capture line. 



2.3 Characteristics of the solar atmosphere during flares from y-ray line measurements 

2.3.1 Atmospheric abundances: deduction from and consequences for y-ray lines 

The fluence in any particular de-excitation line depends on the number and energy distri- 
bution of fast particles, and on the abundance of the species that can contribute to it (either 
directly or via spallation reactions). A set of measured fluences in these lines thus includes 
information on both source chemical abundances, and the numbers and energy distribu- 
tions of the exciting particles. If the energy distributions of the fast particles were known, 
we could deduce a set of source relative abundances from the measured fluences. Solar 
photospheric abundance determinations are fundamental to an enormously wide range of 
astrophysical topics. Determinations in the y-ray source region would complement these, 
extending knowledge of the outer solar atmosphere. This could be particularly important 
at a time when 3-D radiative-transf er modeling has pro duced several significant revisions 
of solar heavy element abundances llAsplund et al.ll2005T) . with the caveat that exceptional 
conditions might occur in the flaring atmosphere. 

If, on the other hand, the source chemical abundances were known, we could deduce 
the energy distributions of protons and a particles in principle. This depends on the ex- 
tent to which the various cr oss-sections have different dependences on fast particle energy 
(Toner & MacKinnon 2004). In practice neither of these sets of variables is known with 
much reliability, and both the abundances and the paramete rs of the energy dis tributions 
should be deduced together. To date, the SMM/GRS work of lMurphv et alj Tl99lh remains 
the only published study where source abundances and particle energy distributions are de- 
termined together. Other stu dies have determine d particle distributions using, e.g., coronal 
or photospheric abundances llRamatvetal.lll995l) . further deciding between these possibili- 
ties using goodness-of-fit tests. Trends in li ne fluence ratios across many events have been 
used to make at least qualitative statements Hh are & Murphy|[T995h ; with many events one 
can attempt to disentangle the separate influences of particle distributions and abundances. 
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Target abundances may be determined non-parametrically or as a choice among several al- 
ternative sets. Ion energy distributions are always represented by assumed analytical forms, 
however (e.g., power-law or Bessel function Ki in momentum) whose parameters are de- 
termined as best fits to data. From these various studies, certain general statements may be 
made (see below) - always noting that these might be subject to change via future revisions 
of nuclear cross-sections or refinements of data analysis techniques. We now survey some 
results of these studies. 

Abundances of several species are enhanced compared to standard phot ospheric values. 
As is found using other measures of coronal abundances (e.g., iMeveril 19931) . enhancements 
may be ordered by first ionization potential (FIP) of the species. Low-FIP elements (Mg, 
Si, Fe) generally show enhancements by an amount that varies from one event to the next 
but may be as great as 3-4. There is a definite "FIP effect," with low-FIP elements show- 
ing more variability than high. At least in SMM/GRS data, the 1.63 MeV line from 20 Ne 
needs both a 20 Ne enhancement (Ne/O = 0.25, rather than the photospheric value of 0.15) 
and a steeply fa lling ion energy distrib ution that can exploit its low threshold for excita- 
tion by protons dShare & Murphvlll995l) . Although an enhanced Ne abundance in the solar 
interior might have h elped with certain problems in understanding helioseismology data 
jDrake & Testa! 2005), there is no convincing evidence for such an enhancement at other 
wavelengths dSchmelz et al. I I200I) . It now seems likely that we are identifying some pe- 
culiarity of the flaring atmosphere rather than the normal state of affairs in the non-flaring 
Sun. 

Rather than a complete fit of the whole spectrum, one often tries to make progress using 
ratios of pairs of line fluences, for instance the ratio mentioned in Section l2.2.5l of fluences 
in the 1.63 and 6.13 MeV lines. Let <P E stand for the measured fluence (cm -2 ) in a line at 
photon energy e (keV). Then the assumed abundances of the emitting species play a crucial 
role. With a photospheric or a coronal abundance ratio for neon and oxygen, one can only 
accommodate the observed values of 4>6.i3, "£1.63 and <Pqa-qs by allowing protons and a s 
to have significantly different energy distributions ( lTonerjf2004l : lToner & MacKi nnonll2004) . 
Because the true degree of enhancement of 20 Ne is uncertain, the proton and a energy dis- 
tributions are correspondingly uncertain. This uncertainty adds to that introduced by the 
different energy-dependence of the cross-sections for excitation of the lines by protons and 
as. For example, assuming a/p = 0.1 among the fast particles, 06.13/^1.63 f° r SOL1981- 
04-27T09:45 (X5.5) implies a power-law energy s pectral index 8 rou ghly in the range 4.3 - 
4.8 if Ne/O = 0.15, but 3.6 - 4.1 if N e/O = 0.25 dRamatv et al.ll 19951) . the former value be- 
lieved appropriate to the photosphere dAsplund et al.l2005l) . Ratios involving the 2.223 MeV 
line give a valuable further constraint in the face of this uncertainty. All target species can 
contribute to free neutron production so some uncertainty in the abundance of any one has 
much less impact on deductions. On the other hand, the intensity of the 2.223 MeV line 
depends on other quantities that must be constrained, if at all, by other means: pitch-angle 
distribution and magnetic geometry in which the fast ions move. 

Usually one needs to integrate over the whole of a flare to get a y-ray spectrum of ade- 
quate statistical quality. There have, however, bee n a few attempts to study the time develop- 
ment of the y-ray spectrum. Mur phy et al] d 19971) found evidence for an increase with time 
of the abundances of low-FIP elements relative to high-FIP ones, from the CGW3/OSSE 
data of SOL1991-06-04T03:37 (X12.0). Evide nce for a similar tr e nd was found in RHESSI 
data for the SOL2002-07-23T00:35 (X4.8) by IShih et all d2003l) . IShare & Murphvl d2006h 
additionally show such a trend in the flare of SOL2003-1 1-03T09:55 (X3.9). It is interesting 
that a similar trend is found in data from three different flares, involving two completely 
different instruments, although the ubiquity or otherwise of this effect remains to be shown. 
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Assuming, unjustifiably, that this effect proves commonplace we may make some prelimi- 
nary comments. We can imagine at least th ree reasons why y-ray source region abundances 
might change (see also iMacKinnonl 120060 . First, abundances might vary with position in 
an active region, and the regions of fast ion precipitation change as the flare progresses. It 
is not at first easy to see why this would produce a systematic change in the ratio of low- 
FIP to high-FIP elements, however. Second, abundances may well vary with height in the 
atmosphere. Even if the atmosphere remains static as the flare progresses, an ion energy dis- 
tribution evolving in time (e.g., hardening) would sample an evolving range of atmospheric 
heights and thus relative abundance s. Third, relativ e abundances might indeed evolve in the 
sou rce region du ring the flare (e.g., IWingleelfl98^ . 

iKiener et al.U2006l) studied the INTEGRAL data from the SOL2003-10-28T1 1 : 10 (X17.2) 
flare, noting a change in time of the relative intensities of the 4.44 MeV 12 C and 6.13 MeV 
16 lines. With line shape data providing additional constraints, their preferred explanation 
involves a change in time of the a/p fast ion abundance ratio, rather than the abundances of 
the target species. This example shows the importance of the extra information obtainable 
from line shapes. 

The time development of the 2.223 MeV line gives a quite different sort of constraint 
on the abundance of 3 He. Indeed, after thermalization (typical duration around 100 s), the 
free neutrons produced in nuclear reactions are captured in the solar photosphere in two 
competitive processes: radiative capture on protons to form deuterium (and the 2.2 MeV 
line) and non -radiative capture on 3 H e ( 3 He(n,p) 3 H). This method was first described and 
employed by IWang & Ramatvl dl974f) . They showed that the greater the 3 He abundance, 
the more rapid the decay of the 2.2 MeV line. The determination of 3 He/ H was there- 
after done on a few flares observ ed by HEAO-1 llHudsonetalJll980l) . SMM dChupp et alj 



198ll: iHua & Lingenfelterill987al) . Gft4iVA7y PHEBUS dTrottet et al.l[l993h . CGRO/OSSE 
Share & Murphyll998h . CGRO/EG RET dDunphv et al.ll999l ; lDunphvll999l). CGffO/COMP TEL 



dRank et al.ll200lh and Yohkoh/GRS dYoshimorilll999al) ; seeTe.g jMandzhavidze & Ramatvl 



2000) for a review of observations before the RHESSI era. A n upper limit on the pho- 
tosphe ric 3 He/H of 1.6 iO.LriO -5 was found in the review of Mandzhavidze & Ramatv 
(2000). The analysis of the time profile of the 2.223 MeV line was also achieved for SOL2003- 
07-23T00:35 (X4.8) observed with RHESSI dMurphv et al]|2003h to get a m easurement of 

He/H. It is one of the factors included in the analysis of iNaiman et alJ (2008). 



2.3.2 Conditions in the positron annihilation region 

Positron production by fl are fast ions has been st udied in detail bv lKozlovskv et all ( fl987l) . 
iMurphv et ai] d 19871) and lKozlovskv et all d2004l) . Positrons may be produced via the decay 
of radioactive daughter nuclei, beta decay of excited states of target nuclei, or via n + pro- 
duction, e.g., in p(p,n7r + X)p'. In the first two cases cross-sections have thresholds mostly 
in the 1-10 MeV energy range, while n + production starts at 200-300 MeV/nucleon. Thus 
positron production, slowing down and annihilation will occur over a wide range of heights 
and the resulting y-rays potentially carry information on deep layers of the atmosphere. 

The very br oad energy range of io ns contributing to positron production is illustrated in 
Figure [Z9l from lMurphv et alj ( 120071) . Effective accelerated ion energies are shown for the 
production of positrons for accelerated 3 He/ 4 He = (white boxes) and 1 (grey boxes) and for 
accelerated a/proton = 0.1 (panel a) and 0.5 (panel b). For very hard spectra dominated by 
protons, the effective ion energies range from 0.1 to > 1 GeV/nucleon because the positrons 
result mainly from very high-energy reactions producing charged pions. At the opposite 
extreme, positron production from soft energy distributions rich in as and 3 He is dominated 
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Fig. 2.9 Effective accelerated ion energies for production of positrons as a function of spectral index for 
accelerated 3 He/ 4 He = (white boxes) and 1 (grey boxes). The effective energy ranges are the 50% yield 
range as defined in Section [2.2.51 Panel (a) is for accelerated a/proton = 0.1 and panel (b) for 0.5. Also 
shown are the effective energies for the 4.438 MeV 12 C line (black boxes) from Figure l2~5l (from Murphy et 
al. 2007). 



by ions of 1-15 MeV/nucleon and even includes a significant contribution from He above 
~1 MeV/nucleon. Note that the effective energies shown in Figure |Z9l are for the production 
of positrons. The escaping 0.5 1 1 MeV annihilation photons can be significantly attenuated 
if the positrons are produced deep in the solar atmosphere. Pion production can be very deep 
since the high energy ions responsible have very long ranges. As a result, even for very hard 
spectra, the annihilation photons that escape can be mainly from decay of radioactive nuclei 
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Fig. 2.10 RHESSI measurements of the 37/27 ratio vs. 511 keV line width for three flares. Time intervals 
for SOL2003-10-28T11:10 (X17.2) are: (a) 11:06:20-11:08:20 UT, (b) 11:08:20-11:10:20 UT, (c) 11:10:20- 
11:16:20 UT, (d) 11:16:20-11:18:20 UT, and (e) 11:18:20-11:30:20 UT. The SOL2003-11-02T17:25 (X8.3) 
measurement was made between 17:16 and 17:26 UT, when the line was broad. The data point when the line 
was narrow is consistent with point (e). The SOL2002-07-23T00:35 (X4.8) measurement was integrated over 
the entire flare. The temperature scale and curves showing t he calculated 3y/ 2y ratio vs. 511 keV line width 
for different densities are for a fully ionized medium (from lShare et all2004l) . 



rather than pions. Such nuclei are produced by ions with energies that are much lower than 
those indicated in the Figure. 

Po sitron annihilation and consequent radiation are dealt with in great detail by Murph y et alj 
( 120051) . Positrons must slow down and thermalize to contribute to 0.5 1 1 MeV radiation. 
They may annihilate directly on free electrons, or via positronium formation. Positronium 
formation can involve both free or atomic electrons, in the latter case via charge exchange re- 
actions. Positronium in turn may exist as para-positronium (net spin 0) or ortho-positronium 
(net spin 1). Para-positronium annihilates to give two photons in the 0.511 keV line (2y 
mode). To conserve spin, ortho-positronium must annihilate producing three photons, at en- 
ergies < 0.511 MeV (3 7 mode). The lifetimes of the two spin states are different, so the 
proportions that annihilate (i.e., the ratio of the two-photon (line) and three-photon (contin- 
uum) spectral components) depend on the rate of collisional destruction and thus on ambient 
density. 

Positrons in an ionized medium of temperature T (K) p roduce a line of FWHM (full 
width at half maximum) width a (keV) llMurphv et alJ2005h 

a = O.Olir 1 / 2 . (2.5) 



For low temperatures an additional, broader component may be discriminated, from positro- 
nium formed via charge exchange in flight. Equation d2.5l l is valid for 0.5 1 1 MeV line emis- 
sion formed via both direct annihilation and para-positronium formation on free electrons, 
but not for positronium formation via charge exchange with atoms or ions. Thus the ioniza- 
tion state of the medium also influences the shape and width of the line. 
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Fig. 2.11 Time histories of the bremsstrahlung, total nuclear de-excitation line, annihilation line and positro- 
nium continuum fluxes, and the width of the annihilation line during SOL2003-10-28T11:10 (X17.2) (from 
IShare et al. 20(3). 



In summary, then, the shape of the line and the relative intensities of the 2y and 3y 
components of the annihilation spectrum include information on the temperature, density 
and ionization state of the positron-annihilation region 

RHESSI has observed 0.511 MeV line emission from the flares SOL2002-07-23T00:35 
(X4.8), SOL2003 -10-28T11:10 (X17.2), SOL2003-11-02T17:2 5 (X8.3) and SOL2005-01- 
20T07:01 (X7.1) dShare etal j2003all2004lMurphv et alj2005h . During the greater part of 
these flares, until quite late times, the width of the line was determined to be 4-8 keV (see 
Figure l2~T0l l. 

The SOL2003-10-28T11:10 flare provided the best-observed 0.511 MeV line. For the 
first 10 minutes of the flare its width is more or less unchanged at ~8 keV. Equation d2.5l > 
then implies an ambient temperature of ~5 x 10 5 K (see Figure [2~TTT >. The measurements of 
the positronium continuum to line ratio (3y/2y ratio) implies an ambient density in the range 
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~ 10 cm (see above) jMurphv et aT1l2005l) . In this initial phase the positron annihi- 
lation region appears characterized b y an unexpected combination of chrom ospheric den- 
sity and transition-region temperature jShare et alj2004lMurphv et alj2005h . T he sim ilarly 
broad line found in SOL2003-10-28T1 1: 10 (X17.2) carries similar implications dShare et aU 
l2004h . The line then narrows in ~2 min, settling down to a width of ~1 keV for the remain- 
ing 12 min of the flare. The best explanation of the late-time narrow line is in terms of 
a medium at a temperature of ~5000 K, but with a relatively high degree of ionization 
(~20%), to suppress the broader feature that results from positronium formation in flight 
(Figure l2~T0l . 

The line in SOL2003-11-02T17:25 (X8.3) is more cons istent with a "normal" solar at- 
mosphere, specifically with the conditions in the layer of the lVernazza et al.l d 198 lh model C 
atmosphere at which the temperature is 5000 K. So far the 0.51 1 MeV line shapes have been 
calculated for single regions with homogeneous conditions. As we see, these yield signif- 
icant insight into at least average annihilation region conditions. There remains a need for 
such calculations to take account of the broad range of depths over which positrons will be 
produced, and the consequent wide range of conditions in which they will annihilate. 

The conditions found for SOL2003-10-28T11:10 (X17.2) are particularly surprising, 
unlike those found at any layer of existing models for quiet-Sun or flaring atmospheres. 
iRavmond et al.l d2007h study lines of O VI UV lines and X-ray emission from SOL2002- 
07-23T00:35 (X4.8). The density and temperature implied by the positron annihilation line 
would have to apply only in a very narrow layer, ~100 m, to be consistent with the measured 
O VI line strengths. Even if such a narrow layer could be produced, bearing in mind the range 
of heights over which fast electrons would deposit energy, we would also have to find some 
mechanism which would localize most positron annihilations there. A multi-wavelength 
study of this flare dSchriiver et al.ll2006h further underlines the difficulty of accounting for 
the positron annihilation line width. 

2.3.3 Inclination of magnetic field lines 

With t he unprecedented spectral resolution in the y-ray domain provided by RHESSI. I Smith et all 
d2003h were able to resolve for the first time the shapes of several of the strongest flare y-ray 
lines in SOL2002-07-23T00:35 (X4.8), determining line widths and Doppler shifts (Fig- 
ure d. 12b (for a complete discussion of the interpretation of line shapes, see Section l2~7ol l. All 
lines were found to be redshifted by amounts that decrease with increasing target species 
mass, as expected for kinematic recoil. The absolute values of these redshifts were larger 
than expected for a flare at S 1 3E72, however, on the assumption of field lines normal to 
the solar surface. ISmifh et al.l d2003h thus suggested that field lines might be inclined to 
the vertical, in this flare at least by ~40°. The alternative explanation, that fast ions are 

highly beamed, appears to be inconsistent with deductions from the a-a line complex in 

J II — k 

the 0.4-0.5 MeV range (Share et al. 2003b), as well as potentially raising difficult questions 
of velocity space stability ( ITamres et aljl 989). Systematic field line tilts to the vertical have 
been suspect e d in the past, on the basis of Zeeman split lines formed at different heights 
dWiehdll978h . iMacKinnon & BrowrJ (H990h had also suggested such a systematic tilt be- 
cause of an apparent asymmetry in the azimuthal distribution of events visible at photon 
energies > 10 MeV. Field lines will fan out with altitude as a result of pressure balance 
( Gabri el 1976h b ut this would produce only a broadening, not a systematic red or blue shift. 

iHarris et all d2007h studied INTEGRAL spectra of SOL2003-1 1-04T19:53 (X17.4) and 
SOL2005-09-07T09:52 (X17.0). These took place on the west and east limbs respectively. 
Systematically tilted loops would produce a systematic difference in Doppler shifts: the 
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Fig. 2.12 RHESSI background-subtracted count spectra for SOL2002-()7-23T00:35 (X4.8). Each panel is 
labeled with the element primarily responsible for the line shown. The carbon and oxygen lines also show the 
secondary peak from the escape of a 5 1 1 keV positron-annihilation photon, which also contains information 
on the line shape. The thick curve shown in each panel is the Gaussian fit with a redshift of respectively 
0.11% for the Fe line, 0.40% for the Mg line, 0.32% for the Ne line, 0.12% for the Si line, 0.79% for the 
C line and 0.58% for the O line plus the underlying bremsstrahlung continuum and broad lines, convolved 
with the instrument response. The thinner li ne is the same fit f orced to zero redshift for comparison. The error 
bars are one o~ from Poisson statistics (from Smith et al. 2003). 



40° tilt suggested bv lSmith et al.l f2003h would imply blueshifted lines from the west limb. 
The two flares studied both display (statistically weak) evidence for redshifts, but none for 
blueshifts. Thus it appears that both flares plus SOL2002-07-23T00:35 (X4.8) involve field 
lines oriented predominantly away from the observer, but that there is no evidence so far for 
a systematic tilt in any preferred sense (Figure [2. 13b . The small number of flares renders this 
study preliminary. At the heart of the interpretation is a picture of unidirectional field lines, 
either normal to the solar surface or all tilted at a single angle to the vertical, instructive but 
very simple. Major y-ray flares originate from complex active regions and will undoubtedly 
involve a broad distribution of magnetic field directions. 



2.4 Narrow lines and ion energy distribution 

Since the first detection of solar y-ray lines in 1972, many solar y-ray line flares have been 
observed with detectors aboard SMM, Hinotori, GRANAT, CGRO and Yohkoh. The analysis 
of these observations have led to quantitative analysis for more than 20 y-ray line events. 
Apart from information on the solar atmospheric elemental abundances which have been 
deduced from the analysis of narrow y-ray lines (see previous sections), energy spectra of 
proto ns have been deduced for 19 SM M/GRS flares ([Share & Murphy Fl995l ; iRamatv et al.l 
1 19951) and for one CGRO/OSSE flare dMurphv et al.(ll997h . Each line is characterized by 



energy-dependent cross-sections for excitation by protons and by a particles. Thus the mea- 
sured fluence in each line provides a sum of energy- weighted measures of the numbers of 
protons and ens (e.g., lRamatvl [l986). above the lesser of the threshold energies for proton 
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Fig. 2.13 Doppler shifts of 12 C and 16 line energies from rest, multiplied by the atomic mass of the recoiling 
nucleus. Dashed line: expected Doppler shift from a simple model of 10 MeV protons incident on ambient 
solar atmosphere nuclei in a downward isotropic beam at 90°. Dotted lines: expected shifts if the protons are 
incident along magnetic field lines tilted at 50° to the horizontal. If this is the explanation for the RHESSI 
point at cos 9 = 0.29 towards the east limb, then the Doppler shifts should follow the branches labelled 
"E" and "W" to the east and west of the Sun's center respectively. Sources of data points are from SMM 
(19 measurements in bins of heliocentric angle; Share et al., 2002), RHESSI meas urements for the flares 
SOL2002-07-23T00:35, SOL2003-10-28 and SOL2003-1 1-02 IShare & Murph\l2006l) . and INTEGRAJJSPl 
measurem ents for the flares S OL20 03-10-28 iKiener et al.l2006l) . SOL2003-1 1-04T19:53 and SOL2005-09- 
07T11:10 lHarris et al.l2007l) : from lHarris et al.l2007l) . 



and a particle excitation. The measured fiuences in more than one line contain information 
on the energy distributions of protons and as. The 1. 63 MeV cross-section o f 20 Ne has a 
lower cross-section than the other de-excitation lines l lShare & Murphvll 19951) , so the ratio 
of the fiuences in the lines at 1.63 and 6.13 MeV, for instance, gives a measure of fast ion 
spectral hardness, representing a variable range of ion energies of around 1-20 MeV (see 
Section [2. 2. 5I >. Higher ion energies, on average, are needed to liberate neutrons in nuclear 
collisions so the ratio of the 2.223 MeV line to one or more of the de-excitation lines sim- 
ilarly gives a measure of spectral shape. The authors listed above assumed the power-law 
form for ion energy distributions and deduced the values of energy power-law spectral index 
from either or both of these measures, for the 20 flares observed by SMM or CGRO. The 
measured ratios imply that the accelerated ion proton spectra should extend as unbroken 
power laws down to at least about 2 MeV/nucleon if a reasonable ambient Ne/O abun- 
dance r atio is used, i.e., in a greement with measurements of the Ne/O abundance ratio in the 
corona dRamatv et all 19951) . The spectrum deduced is thereafter used to estimate the energy 
contained in accelerated protons above 1 MeV. This last determination is, however, largely 
dependent on the value of the ambient Ne/O abundance ratio and also on the composition of 
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Fig. 2.14 RHESSI y-ray count spectrum for SOL200 2-07-23T00:35 (X4.8). The lines show the different 
components of the model used to fit the spectrum (from Lin et al. 2003). 



accelerate d particles, in particular the ratio of energetic He nuclei to energetic protons (a/p) 
(see lRamatv et al.ll 19950 . 

Four y-ray line events so far observed with RHESSI provid e new opportunities for such 
analyses: SOL2002-07-2 3T00:35 (X4.8) (see Figure [27T^T > (s eelLin et alfcOOllSmith et all 
120031 : IShare et alJl2003ah . SOL2003-10-28T11:10 (X17.2) dSfiare et alj|2004T): SOL 2003- 
11-02T17:25 (X8.3) and SOL2005-01-20T07:01 (X7.1) (see lShare & Murphy||200rj , for a 
review). RHESSI has for the first time the energy resolution necessary to resolve all the 
y-ray lines, except the intrinsically narrow 2.223 MeV line, and to determine the detailed 
line shapes expected from Doppl er shifts, thus allo wing deduction of velocity distributions 
of the interacting energetic ions JSmifh et alj|2003l) . The comparison of fluxes in the 20 Ne 
line at 1.63 MeV and of the 12 C and 16 lines at 4.44 and 6.13 MeV was performed for 
these four flares to provide information on the proton spectrum. Apar t from SOL2002 -07- 
23T00:35, for which the proton spectral slope was found to be ~3.5 (Lin et al. 2003 ), th e 
other events have been reported to have much harder slopes (~2.5) |Sh are & Murphvll2 006). 
harder than the average (~4.3) measured previously for the 19 SMM flares. New analyses of 
the most recent events are currently being un dertaken using the improved nuclear continuum 



modeling developed in Murph y et al.il 2009h . which could lead to softer proton spectra than 
s estimated in S hare & Murphvl d2006h . 

The re are also well-reso lved INTEGRAL spectra for three flares: SOL2003-10-28T1 1 : 10 
(X17.2 ( lKieneretal.i r2006). SOL2003-1 1-04T19:53 (X17.4) and SOL2005-09-07T09:52 



(X17.0) jHarris et al 



J I2007I) . The proton spectral index was deduced for these events from 
the ratio of the 2.2 MeV line to the sum of the 4.44 and 6.13 MeV lines. Values of the 
power-law index in the range [3.5, -4] were found. For SOL2003-10-28T1 1: 10 (X17.2), ob- 
served both by RHESSI and INTEGRAL, the difference in the deduced power-law indices 
may result from the different methods used to infer the spectral index f rom observations. 
New analysis using improved nuclear line modeling dMurphv et al]|2009l) could help to un- 
derstand the apparent existing discrepancy. 
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2.5 Alpha particles 

Both protons and as excite most of the strong de-excitation lines. Definitive information on 
the relative numbers of protons and as must be obtained by determining fluences of lines 
to which only one or the other species contri bute. Information on th e a / p ratio had been 
deduced for five flares observed with SMM are & Murphvl 1 1 9981 : iMandzhavidze et ail 
Il999h from the fluence ratio of two lines: the prompt Fe line at 0.847 MeV which is pro- 
duced by the interaction of accelerated protons and as on ambient iron, and a pure a line at 
0.339 MeV which only results from the inter action of energetic a particle s on iron, produc- 
ing an excited state of nickel. It was found by IMandzhavidze et alj ( fl999h that for these five 
flares a/p exceeds the standard value of 0.1 and can even reach 0.5. Nuclear interactions of 
accelerated 3 He with ambient 16 result in three y-ray lines at 0.937, 1.04 and 1.08 MeV. 
While the fluence of the line at 0.937 MeV could be determined for a few flares observed 
by SMM/GRS and CGRO/OSSE, the other two lines cannot be separated from a lines at 
1.05 and 1.00 MeV ( Share & Murphvl [l998h . These lines combine to yield to an unresolved 
feature centered at 1 .02 MeV. Using the information from this unresolved feature and from 
the line at 0.937 MeV, an estimate could be obtained of the ratio of accelerated 3 He with 
respect to accelerated 4 H e. In 7 flares this led to an enhancement of this ratio (0.1 to 1) with 
respect to coronal values ( Mandzhavidze ~et~aT]| 1 999b . In these analyses it was assumed that 
all a ccelerated species have the same energy distributions. As will be discussed in Section 
12771 llbner & MacKinnon! d2004h attempted to show how line fluences might be analyzed 
without this assumption. Finally, observations with good energy resolution help to reach 
stronger conclusions on these issues, and the additional information from y-ray line shapes 
also contributes, as discussed in the next section, to the determination of the a/p ratio. 



2.6 Angular distributions of ions and a/p ratio as deduced from narrow line shapes 

Profiles of narrow y-ray lines can be used to derive the angular distributions of the emitting 
ions (in the interaction site) and of the accelerated ions, if the emission mechanisms are 
coupled with transport models. The nuclear excited states produced by collisions between 
accelerated and ambient particles have indeed very short lifetimes so that the decay y-ray 
is emitted with a Doppler shift due to the recoil from the original collision. Furthermore, 
the collisions with the heavier a particles will produce greater recoil and should produce 
a broad "red" tail to the line shape. Fitting this tail, i.e., looking at the line shape, can 
constrain the a/p ratio independently of line fluxes and ratios. Accelerated nuclei heav- 
ier than He also produce flare y-ray lines but they are largely broadened and merge into 
a quasi-continuum y-radiation see Section [2. lit . The shape of the narrow y-ray lines ul- 
timately depends on the angular distribution and spectrum in the int eraction site, the as- 
sumed a/p ratio, and th e viewing angle of the observer (see, e .g., Ramatv & Crannell 
1976: iRamatv et al.lll979l: iMurphv et alJll988l: IWerntz et al]|l99d : lLang & Werntzlll99ll : 



Kiener et al.ll200ll : lMurphv et al.ll20o'7E 



Before the launch of RHESSI and INTE GRAL, line redshifts and widths had been stud- 
ied for 19 flares using data from SMM dShare & Murphyj [T995l : IShare et alj2002h . Measure- 
ments of energies and widths of the strong flare y-ray lines of C, O and Ne were performed 
as a function of heliocentric angle using data with moderate spectral resolution provided by 
SMM. Flares were grouped in five bins in heliocentric angle, with analysis performed on 
these sums of flares located in the same heliocentric interval (see Figure l2~T3l . Redshifts 
of the Ne, C and O lines of the order of 1% are found for flares close to the disk center 
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Fig. 2.15 Observed and calculated line shapes of the 4.44 and 6.13 MeV ambient I2 C and 16 de-excitation 
lines observed by INTEGRALISPl for SOL2003-10-28T1 1:10 (X17.2). The dashed line represents the calcu- 
lated line shape for a downward isotropic distribution. The full and dott ed lines represent the calculated line 
shapes for the distribution com puted with pitch-an gle scattering as in Murp hy et all 120071) for A=30 (full 
line), A=300 (dotted line) (from lKiener et all2006l) . 

but no redshifts are found for flares close to the limb. The mean widths of the de-excitation 
lines were of the order of 3% and did not exhibit any significant variation with heliocentric 
longitude. It should be noted that the line shapes result from the y-ray spectral analysis and 
are thus not completely independent of the way the analysis is performed. The observed 
redshifts are compared with predictions of line shapes in the case of angular distributions in 
the interacting sites being represented by a downward beam, a fan beam, and a downward 
or upward isotro pic distribut i on. As suming a mean spectral index for ions of ~4.2, a/p=0.3 
and Ne/O=0.25. lshareetal.l J2002h concluded that the measured redshifts as a function of 
heliocentric angle are not consistent with beams in the interacting site but are rather consis- 
tent with a broad angular distribution in the downward direction. These angular distributions 
are of course descriptive of the interaction sites and should be related to the angular distri- 
butio ns in the acceleration region itself by developing transport models (see Murphy et al. 
l2007l) . 

Line widths were clearly ob served for the first time in the flare SO L2002-07-23T00:35 
(X4.8) observed with RHESSI llSmith et alJl2003l : IShare et al.ll2003ah (see Figure l2~T2l . 
It was found that redshifts varied as expected with the mass of the target species, but 
were larger than expected for line production by vertically precipitating ions in a flare at 
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Fig. 2.16 Contour lines in the a/p -A© (width of the injected particle distribution) parameter space of 
downward directed particle distributions for simultaneous fits to the 4.4 and 6.1 MeV lines for different 
phases of SOL2003-10-28T11:10 (X17.2). Confidence levels at 50%, 70%, 90% are plotted for fits using 
an abundance of C reduced by 50% with respect to O as compared to the chromospheric abundance and 
for power-law particle spectra of —3.0 (dotted lines) and —4.0 (continuous lines). The best values of a/p 
and A & for the tw o spectral indices are indicated by open (s = —3.0) and filled symbols (s = —4.0) (from 
iKiener et all2006h . 



this location. The large shifts obtained rather suggest that the magnetic loops along which 
ions propagate are tilted with respect to the vertical (see Section I2.3l l. Two other y-ray 
line flares have been observed with RHESSI for which line shapes have bee n analyzed 
(SOL2003-10-28T11:1 (X17.2) and SOL2003-11-02T17:25 (X8.3) (see, e.g.. iGar] 120(53 ; 
IShare& Murphy] 120061) . The SOL2003-10-28T11:10 (X17.2) event was also observed at 
high s pectral resolution with the INTEGRAL/SP1 spectrometer dGros et alJ200 4; Kiener "et~aT] 
2006), which provided high-resolution data for the 4.44 MeV 12 C and 6.13 MeV 16 lines 
(see Figure l2~15l l. The comparison of the results of detaile d calculations of lin e shapes to the 
observations places strong constraints on the a/p ratio dKiener et alj|2006h . In particular, 
the observed line width is inconsistent with an interpretation where too many a particles 
relative to the number of protons would be produced or otherwise the lines would be too 
broad when compared with the observations (Figure [2.15b . However, the line shapes depend 
on many parameters: the angular distribution of the interacting ions, the spectral index of the 
energetic ions, and the a/p ratio; these cannot be determined independently so the analysis 
results in extended regions of allowed parameter space. 

Using a spectral inde x between 3 and 4 de duced from the ratio of the 2.2 MeV to the flu- 
ence at 4.4 and 6.1 MeV dTatischeff et alj2005h . simultaneous fits of the C and O line shapes 
and line fiuences have led IKiener et al.l ([2006) to favor the production of the y-ray line emis- 
sions in the SOL2003-10-28T11:10 (X17.2) flare, by downward directed ion beams with a 
relatively low a/p ratio, around 0.1 (Figure l2.16l l. High-resolution y-ray line spectros copy 
of two other weaker flares located at the limb was also performed by Harris e t al. (2007j) us- 
ing INTEGRAL data (SOL2003-1 1-04T19:53 and SOL2005-09-07T09:52; see also Section 

In addition to the line width of the narrow flare y-ray lines discussed above, the profile 
of a/ 4 He fusion lines around 0.452 MeV (a/ a lines) may also provide a potential diag- 
nostic of the angular distributions of energetic particles in the interaction site. From SMM 
observations of two flares, it had been shown that the line shape of the a/ 4 He fusion line 
was inconsistent with its production by a downward beam and was rather consistent with 
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Fig. 2.17 Count spectrum of the a- 4 He line measured by RHESSI for the SOL2002-07-23T00:35 (X4.8) 
event after subtraction of the other spectral components. The curve shows the RHESSI response to the line 
shape in the case of an a particle distribution with saturated pitch-angle scattering (see Section l2.2.U prop- 
agating in magnetic loops perpendicular to the solar surface at a heliocentric angle of 73° (from Shar e et al l 
2003b). 



an isotropic or a fan-beam distribution (i.e., a s obtained at a magnetic mirror point) in the 
interacting site (Share &Murphy|[T99lll998l) . This is consistent with what was found by 
IShare et aljfcooii ' by analyzing line shapes for SMM flares. 

Since the launch of RHESSI, high-resolution spectra of the a/a lines have been obtained 
for the SOL2002-07-23T00:35 (X4.8) (Figure ITTTI l and SOL2003-10-28T11:10 (X17.2) 
events dShare et al] |2003d . l2004h . For these two flares the shape of the line complex is also 
found to be consistent with a downward isotropic distribution. 



2.7 Energy content in ions 

De-excitation y-ray lines provide information about ions above ~1 MeV/nucleon. The en- 
ergetic ion distribution below this energy is essentially unknown but is of particular inter- 
est for the total ion energy content (see Section 12. 10b . Since the first detection of y-ray 
lines in 1972, quantitative results from y-ray spectroscopy have been obtained for more than 
20 events, yielding information on the energy content in ions. Several ratios of pairs of nu- 
clear lines have been used to i nfer the ion energy s pectrum and then the ion energy content: 
&2. 223/^4-7 and ^1.63/^6.13 teamatv et al.| [l996). As an example, the threshold energies 
for excitation of Ne and O lines by fast protons are significantly different, >2 and >8 MeV, 
respectively, which allows the use of the ratio of their measured fluxes to infer the ion energy 
distribution ( Share & Murphv 1995). This determination is not independent of the assumed 
target abundances but combining results from several such methods reduces some of this 
uncertainty. The comparison of the steepness of the ion spectra deduced from the different 
ratios, ^2,223/^4-7 and ^1.63/^6.13 f° r different flares, shows indeed some consistency only 
when an enhan ced Ne abundance is assumed, Ne/O = 0.25 instead of the more standard 
value of 0.14 dRamatv et al.lll996h (see also Section [2.3.1t . As most of the narrow lines can 
be produced by energetic protons and as, the usual hypothesis is that the energy distributions 
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of accelerated protons and as expressed in energy/nucleon are identical. The relaxation of 
this assumption will be discussed below. The determination of the proton spectrum from the 
ratio of the excitation lines from Ne and O depends on the assumption made on the ratio of 
Ne/O in the flarin g region and on the ratio a/p of fast a particles to energetic protons. It 
is indeed found in lRamatv et al.l Jl995h that steeper ion spectra are needed to reproduce the 
same ratio of line fiuences for smaller values of Ne/O and larger values of the a/p ratio. The 
analysis based on the 19 events observed by SMM 7GRS led to the conclusion that the energy 
contained in > 1 MeV ions lies in the ranee 

10 29-32 erg for 

y-ray line flares (see Ramaty et 
al. 1995, 1996, and the summary by Miller et al. 1997). These estimations were done un- 
der the following assumptions: the slope of the ion power-law spectrum is deduced from 
the ratio of the Ne/O line assuming thick-target emission and values of Ne/O=0.25, a/p 
=0.1 and "impulsive" composition for other species (i.e. increase of Ne/O, Mg/O, Si/O, and 
S/O by a factor of ~3 and by a factor of ~10 for Fe/O with respect to coronal composition 
for the accelerated particles). The power-law spectrum is extended down to 1 MeV with 
a flat extension to lower energies. For the giant flares SOL1991-06-01T16:14 (X12) and 
SOL1991-06-04T03:37 (X12.0), observed either by GRANAT/PUEBVS and CGRO/OSSE, 
the energy contained in 1 MeV/nucleon ions has also been de duced and is found to be in the 
10 32 -10 33 ergs range llRamatv et ail 19971 : iMurphv et ail 1997b . ForSOL1991-06-04T03:37, 
the total energy contained in accelerated protons is 1.7 x 10 32 erg and the total energy con- 
tained in accelerated ions is 1.0 x 10 33 erg under the following assumptions: impulsive flare 
abundances for accelerated particles, identical energy dependence for all ion spectra, and 
a/p = 0.5. 

The relaxation of the assumption of similar p roton and a energy spectra in MeV/nucleon 
has been studied bv lToner & MacKinnon! j2004l) and applied to the estimates of the ion en- 
ergy content in flares. As some of the narrow lines (in particular the a-a lines at 0.429 
and 0.478 MeV) result only from the interaction of energetic as with the solar atmo- 
sphere, the measured fiuences of these lines can be used to infer the number of a particles 
>1 MeV/nucleon once the steepness (8 a ) of the a spectrum has been chosen. This allows a 
computation of the flux of the narrow lines at 1.63 MeV and 6.13 MeV produced by the as. 
The remaining fiuences in the narrow lines should be produced by fast protons. The ratio 
of the remaining fiuences is then used to determine the steepness (8 P ) of the proton spec- 
trum and then the number of fast protons required to produce the remaining fiuences. Note 
that some values of the steepness of the a spectra (excessively hard ones) can be excluded 
directly by the fact that the fiuences of the produced narrow lines at 1.63 and 6.13 MeV 
would exceed the observed ones. The reexamination of the spectra and numbers of parti- 
cles d educed from SMM/GRS spectra of four flares chosen in the list of IShare & Murphv 
dl995h shows that harder spectra for a particles than for protons are generally deduced, that 
most of the energy still goes to protons (by a factor of 100 to 1000) and that as a whole 
less energy goes to ions. As an example, for the well studied SOL1981-04-27T09:45 (X5.5) 
event, a total i on energy of 1.6 x 10 31 erg is found with 8 a =2.5, 8 P =4A, a/p=0.004 from 
the method of iToner & MacKinnon! d2004b . while the usual analysis (using the ratio of the 
fiuences of the 1.63 MeV and 6.13 MeV lines) would lead, for the usually chosen value of 
a/p= 0.5, to a total ion energy of 5.5 x 10 31 erg with an a and proton spectral index of 5.2. 
This discussion clearly shows the present limits and the uncertainties in deducing ion energy 
c ontent from y-ray line spectro scopy. It must finally be noted that in the analysis developed 
in Toner & MacKinnon (2004) the meaning of a/p is different from the standard one when 
both ion spectra are similar. In the case of different ion spectra, the value of a/p which is 
indicated above is the total number of a particles > 1 MeV/nucleon to the total number of 
protons above 1 MeV. 
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Other effects can also alter the amount of energy estimated for accelerated ions as de- 
duced from y-ray line spectroscopy. As stated above, the production of the y-ray line emis- 
sion is usually estimated under the assumption of thick-target production. This means that 
the energy loss rate of the radiating ions is also a parameter (even if well hidden) of the y-ray 
line production. Most of the models of y-ray line production in flares assume that the target 
in which the fast ions interact and produce lines is a "cold" neutral target (i.e., energetic par- 
ticle speeds much larger than the target particle speeds). This may no longer be true if ions 
with energies close to the threshold energies for y-ray line production interact in a hot flaring 
plasma (temperature of a few 10 7 K), i.e., if the y-ray line emission site is partl y coronal (see 
Sectio n |5]for some discus sion on possible coronal sites of y-ray line emission). lEmslie et al.l 
dl997h . and more recently iMacKinnon & Toned y003), have examined the effects of intro- 
ducing a "warm target" for the production of y-ray lines. In such a warm medium, fast ions 
with energies just above the threshold for nuclear line production will lose less energy to 
collisions than in a cold medium. This has the effect of increasing the efficiency of line pro- 
duction by ions just above the threshold. This increase is counterbalanced by the fact that 
the energy losses in an ioni zed (because hot) medium are larger than in a neutral medium. 
IMacKinnon & Tonerj d2003l) estimated the strength of the Ne 1.63 MeV line relative to the 
strength of the O 6.13 MeV line, used to estimate ion spectra, as a function of the temper- 
ature of the target and as as function of the a/p ratio, i.e. assuming as in most works a 
similar spectral shape for energetic protons and as in MeV/nucleon. They show that in a 
warm target, there is an increase of the relative production of the Ne 1.63 MeV line to the 
O 6. 13 MeV line. This increase is more important for steeper spectra and for lower values of 
a/p (this simply illustrates the fact that the Ne line cross-sections have much lower thresh- 
olds for proton production than for a production). Flatter ion spectra and thus lower ion 
energy contents are then required to reproduce the observed line fluences ratio than usually 
deduced (a reduction by more than a factor of ~10 of the total energy contained in protons 
and as can be achieved if the emission site has a temperature of 6 x 10 7 K and a/p is 0.1; 
see MacKinnon & Toner, 2003). Taking into account the effects of a warm target could fi- 
nally reduce the discrepancies which were found in the ion spectral indices deduced from 
the different line ratio ^2.223/^4-7 and <P\ 63/^6.13 f° r a standard value of Ne/O. Indeed, 
even with a standard value of the Ne/0 abundance ratio, the inc lusion of a warmer emissio n 
site leads to harder ion spectra to produce the same ^1.63/^6.13 ( MacKi nnon & ToneJ2003l) . 



Observations performed by RHESSI brought information on the ion energy content for 
three more y-ray line events: SO L2002-07-23T00:35 (X4.8) dLin et alj|200a: lEmslie et al.1 
2004d : lDauphin & Vilmerl2007l) . SOL200 3-10-28T1 1 :10 (X17.2) dGarJ2005| : lshare & M urphy 
20061) . and SOL2003-1 1-02T17:25 (X8.3) jGanl20bl : IShare & Murprrvhooa) . Usingapro- 
ce dure similar to the on e wh ich has been used pr e viously fo r SMM/GRS and CGRO events 
bv lMurphv etafl J 19971) and lRamatv et all dl997l) . lLm et al.l J2003I) found a minimum total 
energy in accelerated protons above 2.5 MeV of 1.4 x 10 30 ergs (10 31 erg for all ions) for 
SOL2002-07-23T00:35 (X4.8) (assuming a/>=0.5). Consistent values of the total energy 
in protons and a particles above 2.5 MeV/nucleon were found bv lDauphin & Vilmerl ( 120071) 
from the modeling of the y-ray line time profiles. Values of the proton energy content derived 
for S OL2q 03-10-2 8Tll:10 (X17.2) are con sistently found to be around 1.5 x 10 31 ergs by 
iGar] (2005) and bv lShare^ Murphv (2006). For SOL2003-11-02T17:25 (X8.3) the proton 
energy content ranges from 0.5 x 10 31 ergs to 4.2 x 10 31 erg depending on the value of the 
proton spectral index used (3 — 3.6). This clearly shows the still-large uncertainties in the 
determination of these quantities. 
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2.8 Ion and electron acceleration in flares and variability 

Using the accelerated ion and electron spectra derived from y-ray and X-ray spectroscopy, it 
is possible to comp are the energies contained i n flare -accelerated ions and electrons. Such a 
study performed by Ramatv & Mandzha vidzel d2000l) for the 19 y-ray line flares observed by 
SMM showed that the energy contained in > 1 Me V ions may be comp arable to the energy 
contained in subrelativistic electrons (>20 keV). Emslie et a j] (l2004ah obtained a similar 
result for the first y-ray line flare observed by RHESSI. At present, such comparisons of 
subrelativistic electron and ion energy contents have not been performed for more flares. 
The published comparisons show that there is a large dispersion of the relative electron 
(>20 keV) and ion (>1 MeV/nucleon) energy contents from one flare to the other, the 
energy contained in the ions sometimes exceeding the energy in the electrons. However, it 
must be recalled that all these estimations strongly depend on the low energy cut-off of both 
energetic ele ctrons and ions whi ch is for electrons not always easy to determine (see the 
discussion in Holm an et alj|201 lh . 

Another way of dealing with the question of the relative ion and electron acceleration in 
solar flares is to investigate the relationship between ion acceleration diagnosed by y-ray line 
fluences and relativistic electron acceleration diagnosed by the >300 keV bremsstrahlung 
radiation. A first comparison was carried out on SMM/GRS data between the >300 keV 
bremsstrahlung radiation and the excess rate in the 4-8 MeV range (i.e., the rate above the 
bremsstrahlung continuum which is attributed to the production of y-ray lines by ions (see 
Figure 12. lb ) showing that there is a good correlation between ion productio n (normalized 
above 30 M eV/nucleon) and electron production above 300 keV (see e.g. lForrestlll983l : 
Clump 1 19841 ). The relationship between energetic ion and relativistic electron acceleration 
can also be studied by comparing the 2.223 MeV line fluence to the >300 keV X-ray emis- 
sion (see, e.g., Murphy et al., 1993). Since the launch of RHESSI, 29 flares have been ob- 
served with significant > 300 keV continuum fluence. The fluences of the 2.223 MeV lines 
have been estimated for these flares from spectral analysis and corrected for attenuation for 
flares close to the limb. The neutron-capture line fluences are found to be highly correlated 
to the >300 keV bremsstrahlung fluences for all the flares located at heliocentric angles 
less than 80° and over a range of more than three orders of magnitude in fluence. A simi- 
lar correlation was found for the flares previously observed by SMM. The ion energy range 
which produces the 2.223 MeV line depends in fact on the ion spectral index and on the 
accelerated and ambient abundances (see Section [233] and Murphy et al., 2007). As a first 
approximation, this line can be considered as a good indicator of the number of >30 MeV 
protons and its fluence is taken as being proportional to the number of protons above this 
energy. When both SMM and RHESSI data sets are taken into account, the strong correla- 
tion obtained between the neutron-capture line fluence (corrected for attenuation) and the 
bremsstrahlung fluence (>300 keV) (see Figure l2~T8t clearly shows a link between the total 
ener gy content in proto ns above 30 MeV and the total energy content in electrons above 300 
kev lshihetai]2 009b). This suggests that high-energy electrons and ions are directly linked 
by acceleration processes. This may be a more direct link than for subrelativistic electrons 
and ions, given the larger variation between electron energy contents above 20 keV and ion 
energy contents above 1 MeV/nucleon discussed above. This is a constraint which should 
be borne in mind when discus sing accele r ation mechanisms in solar flares. 

The correlation found by IShih et all d2009bl) between ion and relativistic electron ac- 
celeration seems to indicate that energetic ions are produced as soon as there is a signif- 
icant production of energetic electrons above, e.g., 300 keV. From this correlation, there 
does not appear to be a distinct class of electron-dominated flares that produce much more 
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Fig. 2.18 Neutron-capture line fluence as a function of the electron bremsstrahlung fluence above 300 keV 
for flares with heliocentric angle <80° (RHESSI in solid symbols, SMM in open symbols). Circles (triangles) 
represent flares with complete (incomplete) coverage. Colors represent the GOES soft X-ray (SXR) duration 
in three bins. The dotted line indicates the best fit in linear space that passes through the origin with a ratio 
of 0.06 6, and the dashed lines have slopes that differ by factors of ~2 from the best-fit line (from Shi h et all 
[2009bt) . 



high-energy bremsstrahlung radiation relative to nuclear emission (e.g.. lRiegeret 1d Il998h 
although some episodes of the flares reported bv lShihetalJ j2009b) may exhibit higher rel- 
ative numbers of electrons with respect to ions than the mean value over several flares as 
well as over the whole event. Extreme cases of the variability of the electron bremsstrahlung 
component to the nuclear line component have been reported in the SMM era during short 
duration (a few seconds to a few tens of seconds) transient bursts observed above 10 MeV at 
any time during a flare. They are referred to as electron-dominated events (Rieger, 1994) and 
are characterized by weak or no detectable y-ray line emission and by hard > 1 MeV elec- 
tron bremsstrahlung spectra. They were first reported from SMM/GRS obs ervations (e.g.. 



Rieg er & Marschhauserl|l99lf) . and were afterwards observed by GAMMA-1 dAkimov et alJ 
1992|: lLeikov et al.| l99J) ."pHEBUS an d SIGMA experimen ts aboard GR ANAT dPelaez et all 
1992l ; lvilmer et al.lll994l) , CGRO (e.g jDingus et al.lll994l) and Yohkoh ( lYoshimorill 1999b 



The s pectral analysis of 12 electron-dominated events observed by SMM/GRS dRieger et alJ 
Il998h confirmed the hardness of the bremsstrahlung spectra above 1 MeV (mean value of the 
power-slope around —1.84). The mean value of the spectral index between 0.3 and 1 MeV 
(~2.7) does not differ significantly from that of other flares. However, the apparent lack 
of y-ray line emission does not rule out a sim ultaneous production of relativistic electrons 
and ions llTrottet et al.lll998l : lvilmer et al.ll 19991) for the electron-dominated events observed 
G^?A/VA77PHEBUS. Indeed, if one assumes that the energy content in ions above 1 MeV 
is similar to tha t conta ined in electrons ( i .e., aro und a few 10 29 ergs for the cases studied 
by iTrottet et all i ll 9981) and IVilmer et all i 19991) . as appears to be the case for y-ray line 
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flares, one finds that no detectable y-ray line fluence could be observed, given the hardness 
of the electron bremsstrahlung component and the limited spectral resolution (and thus line 
sensitivity) of the PHEBUS experiment. The only remaining question would then be to un- 
derstand why the spectrum of accelerated electrons is so much harder in these events than at 
most other times. 

Variability of the accelerated electron-to-proton ratio not only occurs from flare to flare 
as discussed above , but also on time scales of tens of seconds within individual events. 
IChu pp et al.1 dl993h reported, e.g., for an SMM event a variation of the relative number of 
electrons above 500 keV and of ions above 30 MeV on time scales of tens of seconds. 
iDauphin & Vilmerl J2007h also reported variation of the relative number of electrons above 
150 keV and of the number of energetic protons and a particles above 1 MeV/nucleon in 
the evolution of SOL2002-07-23T00:35 (X4.8) observed by RHESSI. These variations were 
deduced from the analysis and modeling of time delays between y-ray line emissions and 
hard X-rays. A question of obvious interest would be to determine whether these changes 
are connected to spatial source evolution. Such variations of the ratio of accelerated elec- 
trons to accelerated protons should also be taken into account when discussing acceleration 
processes. 

Another interesting result follows from the statistical analysis of lShih et al. (2009b). For 
eight of the flares with the largest 2.223 MeV line fluence, a correlation is found between the 
2.223 MeV line fluence (i.e., the number of >30 MeV protons) and the fluence above 50 keV 
as well as with the GOES class. However, this correlation is found only when a threshold in 
the production of ions is re ached. A similar relationship had been discussed at the time of 
SMM bv lCliver et al. (1994). This ultimately suggests that while the acceleration of protons 
above 30 MeV is closely related to the acceleration of relativistic electrons, the acceleration 
of subrelativistic electrons is only proportional to the acceleration of relativistic electrons 
and ions when a given threshold of high energy particles is reached. This may suggest two 
acceleration processes, one producing proportional quantities of relativistic electrons and 
ions and the other one producing mostly subrelativistic electrons. 



2.9 Long-lived radioactivity from solar flares 

Interactions of flare-accelerated ions with the solar atmosphere can synthesize radioactive 
nuclei, whose decay can produce observable, delayed y-ray lines in the aftermath of large 
flares. The detection of these delayed lines would provide new insights into the spectrum 
and composition of the flare fast ions. Many of these radioisotopes would be produced pre- 
dominantly by interactions of fast heavy ions with ambient hydrogen and helium, especially 
since accelerated heavy nuclei are believed to be signi ficantly enhanced compared to the 
ambient medium composition (e.g., iMurphv etalll99ll) . 

One of the most promising of such lines is at 0.847 MeV resulting from the first excited 
state of 56 Fe which in turn has been produced via the decay of 56 Co (half-life T\ / 2 =77.2 days) 
llRamatv & Mandzhavidze1l2000l : [Kozlovskv et alj2002l:lTatischeff et alj2006u . Another line 



produced in this decay is at 1.238 MeV dShih et alj [2007). For a few days after a large flare 
however, the most intense delayed line is predicted to be the 511 keV positro n-electron anni- 
hilatio n line resulting from the decay of several long-lived jS + radioisotopes. iTatischeff et all 
(2006) have pointed out other delayed y-ray lines that appear to be promising for detection, 
e.g., at 1.434 MeV from the radioactivity of both the isomer 52 Mn m (7^=21.1 min) and the 
ground state 52 Mn« (J 1/2 =5.59 days), 1.332 and 1.792 MeV from 60 Cu (7/ 1/2 =23.7 min), 
1.369 and 2.754 MeV from 24 Na (7/ 1/2 =15.0 hours), and 0.931 MeV from 55 Co (Ty 2 = 
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Fig. 2.19 Total fluences of the delayed lines at 0.847 MeV (solid curves) and 1.434 MeV (dashed curves) 
as a function of the accelerated a/p abundance ratio (see text), for spectral indices of the accelerated ion 
power-law spectrum s=2, 3.5 and 5. The calculations are normalized to a total fluence of 300 photons cm -2 
emitted during the y-ray flare in the sum of the 4.44 and 6.13 MeV ambient 12 C and ls O de-excitation lines 
Cfrom lTatischeff et all2006l) . 

17.5 hours). The delayed lines should be very narrow, because the radioactive nuclei are 
stopped by energy losses in the solar atmosphere before they decay. Unless the flare is very 
close to the solar limb, the y -ray line attenuation in the solar atmosphere should not be sig- 
nificant fsee lHuaetaI1ll989h . as long as the radioactive nuclei do not plunge deep into the 
solar convection zone. Sear ches for upper l imits of the radioactivity line at 1.238 MeV were 
performed with RHESSI ( Shih et al. 2007) for periods after the large flares of July 2002, 
October-November 2003 and January 2005. No strong conclusion has been drawn so far. 

In addition to y-ray lines emitted from de-excitation of daughter nuclei, radioactive X- 
ray lines can be produced from the decay of proton-rich isotopes by orbital electron capture 
or the decay of isomeric nuclear levels by emission of a conv ersion electron. The s trongest 
delayed X-ray line is predicted to be the Co Ka at 6.92 keV dTatischeff et al.ll2~006h . which 
is produced from both the decay of the isomer 58 Co m (7 1/ / 2 =9.04 hours) by the conversion of 
a K-shell electron and the decay of 57 Ni (7^2=35.6 hours) by orbital electron capture. The 
attenuation of this line by photoelectric absorption in the solar atmosphere should be <10% 
for flares occuring at low heliocentric angles, as long as the radioisotopes produced in the 
chromosphere and upper photosphere are not transported to greater depths. Distinguishing 
this atomic line from the thermal X-ray emission can be challenging until the flare plasma 
has significantly cooled down. However, a few hours after the flare end time the thermal 
emission will be gone, or significantly reduced, and the delayed Co Ka line will be more 
easily detected. 

To illustrate the potential of these lines for revealing details of the heavy ion energy dis- 
tribution, we show in Figure l2~19l calculated fluences of the 0.847 MeV and 1.434 MeV lines 
as a function of the accelerated a/p abundance ratio. The calculations were normalized to a 
total fluence of the summed 4.44 and 6.13 MeV prompt narrow lines of 300 photons cm~ 2 , 
which is the approximate fluence observed in the SOL2003-10-28T11:10 (X17.2) flare. We 
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assumed for the abundances of fast C and heavier elements relative to g-pa rticles the aver - 
age composition of solar energetic particles measured in impulsive events llReames| [T999). 
Thus, the predicted fluence variations with accelerated a/p actually show the relative contri- 
butions of reactions induced by fast protons. The fluences decrease for decreasing a/p ratio 
(i.e., increasing proton abundance), because, for a/p > 0.05, the radioisotopes are predom- 
inantly produced by spallation of accelerated Fe, whereas the ambient 12 C and 16 lines 
largely result from fast proton interactions. This effect is less pronounced for s=5, because 
for this very soft spectrum, the contribution of a -particle reactions to the prompt line emis- 
sion is more important. Thus, a concomitant detection of the 0.847 MeV and 1.434 MeV 
lines would allow measurement of not only the abundance of accelerated ions, but also of 
their energy spectrum. 

As well as serving as diagnostics of the flare fast ion population, the radioisot opes can 
serve as tracers to study solar atmospheric mixing dRamatv & MandzhavidzdfeOQOb . A mea- 
surement of the decay profile with time of several delayed lines would place very useful 
constraints on the extent and timescale of mixing processes in the outer convection zone. A 
future observation of the size and development of the radioactive patch on the solar surface 
would furthermore provide unique information on both the transport of flare-accelerated 
particles and dynamics of solar active regions. Note that solar radioactivity can be the only 
way to study flares that had recently occurred over the east limb. 



2.10 Low-energy protons 

The total energy content of a power-law particle energy distribution, F(E) oc E~ s , is domi- 
nated by the lowest energy E for which this form holds (as long as s > 2). The low-energy 
form of the distribution is thus important for assessing the ions' overall importance in flare 
energetics and energy transpo rt. Power-law ion dis tributions measured in space appear un- 
broken to 0.02 MeV/nucleon dReames et al.lll99~7h . If this were to be true also at the flare 
site, extrapolation to these low energies of the ion distributions deduced above ~2 MeV 
dMurphvetalJll997h would yield estimates of ion energy content well in excess of any 
other measure of flare total energy. Unfortunately, there are few diagnostics bearing on ion 
distributions below 1 MeV. 

Proton-capture cross-sections typically have resonances in the 0.1-1 MeV energy range. 
Since these resonances result from formation of the compound nucleus in excited states, 
they also give rise to y-ray lines as these states decay to lower-lying energy levels. The cross- 
sections for these lines are small compared to those for de-excitation, but usef ul constraints 
on to tal ion energy could result from upper limits at the 10~ 5 cm~ 2 level jMacKinnonl 
Il989h . The two strongest lines are at 2.37 MeV, from 12 C(p, 7) 13 N, and at 8.07 MeV, from 
13 C(p,y) 14 N, constraining the proton energy distribution above 0.46 and 0.555 MeV, respec- 
tively. Attempts to constrain the fluxes in these lin es observationall y have proven inconclu - 
sive, however, whether in flares dShare et alj200ll) or the quiet Sun (McConne ll et al J 1997b . 
The line at 8.07 MeV would suffer less competition from the strong de-excitation lines, but 
attempts to constrain it would still require fine energy resolution and/or low instrumental 
background. 

Other possible diagnostics for lower energy ions include Doppler-shifted Lyman-a line 
emission dCanfield & Chandll985h and line impact polarization of Ha, HjS and other lines 
llHenoux et al J 1 1990: Vogt & Henoux H"l996b . Ob servations of Ha and H/3 linear polariza- 
tion in flares seem increasingly well-established dXu et alj|2005l, 120061 : iFirstova et alj|2008b 
although there remains an ambiguity in its interpretation between ion beams and stream- 
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ing electron distributions (e.g., return currents). W hile there is a claime d observa t ion of 
Doppler-shifted Lyman-a emission in a stellar flare dWoodgate et al.lll992h . lBrosiusl feOOlh 
finds no evidence for a-particle beams in He II Lyman-a from a C3.8 flare. 

Their importance for the flare energy budget notwithstanding, flare-site ions below 1 MeV/ 
nucleon remain more or less unconstrained by observations. A promising new technique 
could be to detect suprathermal ions through charge exchange with the ambient hydrogen 
(see the first detection of energetic neutral hydrogen atoms in a solar flare by Mewaldt et al. 
2009). 



2. 1 1 Heavy energetic ions as diagnosed from broad lines 

The relative abundances of accelerated ions heavier than He can be determined from the rela- 
tive fluences of the broad de-excitation lines. These lines are produced when heavy energetic 
ions interact with ambient H and He. Because the heavy ions lose little of their kinetic en- 
ergy in the reactions, the excited nucleus produced has a large velocity and the de-excitation 
y-ray line can be strongly Doppler-broadened by up to ~20% FWHM. Because the lines are 
so broad, many of them overlap, making the heavy-ion abundance determination much more 
difficult than the corresponding determination of ambient abundances using the narrow de- 
excitation lines. In addition, such broad lines are difficult to distinguish from the underlying 
continua formed by the electron bremsstrahlung and the nuclear continuum. The nuclear 
continuum is composed of hundreds of relatively weak, closely spaced de-excitation lines 
that also result from the nuclear interactions along with the strong lines clearly visible in 
flare spectra. These lines are so numerous and closely spaced that they merge and appear as 
a continuum to y-ray detectors. This continuum, however, does have structure, on the order 
of several hundred keV, and knowledge of this structure is critical to reliably separate it from 
the broad lines. The early attempts to derive heavy-ion abundances from flare y-ray data that 
we discuss here used relatively crude modeling of this nuclear continuum. On-going anal- 
yses using improved nuclear-continuum modeling recently obtained from modern nuclea r 
rea ction codes will impr ove the reliability of such determinations (see 
iMurphv et all dl99ll) u sed the y-ray de-excita tion line code from 
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(see the revised version in iKozlovskv et~aT]l2002h to calculate both narrow and broad line 
spectra from the most abundant elements in the solar atmosphere, and the most abundant 
accelerated ions. The abundances were varied to obtain the best fit to the SMM data for 
SOL1981-04-27T09:45 (X5.5). The resulting accelerated-ion abundances relative to ac- 
celerated carbon are shown in Figure 12.201 Also shown for comparison are similar ele- 
ment abundances measured in space from large proton flares (LPF) and impulsive, 3 He-rich 
events. The uncertainties for most of the elements are quite large, but Mg and Fe are seen to 
be significantly (>3 — 4 a) enhanced relative to C and O, similar to the enhancements seen 
in the impulsive 3 He-rich events. 

Obser vations of the beh ind-the-limb flare SOL1991-06-01T16T4 (X12) with GRANATI 
PHEBUS llBaratetal.ll994h have moreover shown that the enhancements in heavy ions may 
increase with time in the course of the flare, r eaching towards its end the highest valu es ob- 
served for solar energetic particles in space dTrottet et a l. 1996; Ramatv et alJI 19971) . As a 
thin-target production of the emission is required to account for the very high observed ratio 
of 1.1-1.8 MeV flux to 4.1-7.6 MeV flux, the temporal evolution of the abundances of accel- 
erated ions is to be related to the evolution of the accelerated particles themselves. It must 
finally be noticed that, although a behind-the-limb flare, this event is associated with one 
of the largest y-ray line fluences observed so far. This flare is also surprisingly associated 
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Fig. 2.20 Derived accelerated-i on abundances rela tive to accelerated carbon (circles) for SOL 198 1-04- 
27T09:45 (X5.5) (adapted from Murphv et al. 1991). Also shown are large proton flare (LPF, squares) and 
impulsive 3 He-rich event (diamonds) abundances (Murphy, private communication). 



with the observation of a strong flux of neutrons by the OSSE experiment dMurphv et alj 
1 19991) .As the neutrons are also expected to be produce d as a thin target, a hard ion spec- 
tral index (power-law index around —2) is deduced by iMurphv et al.l 1 19991) up to at least 
50 MeV. 



3 Pion-decay radiation in solar flares 

Observations of pion-decay radiation from solar flare y-rays and of neutrons (see next 
section) combined with observations of flare y-ray lines give a complete picture of the 
accelerated ion distribution above a few MeV/nucleon to several GeV/nucleon (se e, e.g., 
lRamatvll98&IChupplll984ll99d : IVilmer & MacKinnon1l2003l : lc"hupp & Rvanll2009l. for re- 
views ). Studies of high energy emissions from the flare si te started with SMM (Murphv et al.l 
19871), continued with events observed with GAMMA - 1 jAkimov et al.ll 199% , CGRO (e.g., 
Mandzhavidze & Ramatv lll992l: iKanbach et al.l 1 19 93: Dunphv et al. 1999), and now with 
CORONAS-F dKudela et al.ll2003l) . Pion-decay radiation together with observations of neu- 
trons provide information on the highest ion energies produced during flares. Charged pions 
decay to yield electrons and positrons, which in turn produce y-rays by bremsstrahlung. 
Positrons also contribute to the continuum by annihilating in flight. Neutral-pion (7T°) de- 
cay results in two photons, with one emitted at a high energy. This results in a very flat, 
broad "bump" feature which has a maximum at 67 MeV. Synchrotron losses of electrons 
and positrons may be important and shorten their lifetimes, therefore reducing t heir contri- 
bution with respect to the radiation from neutral-pion decay (Mur phv et all 19871) . Energetic 
electrons above 10 MeV are also produced in solar flares and produce bremsstrahlung con- 
tinuum, potentially masking pion-decay radiation. 

Pion decay radiation was first observed with the SMM/GRS instrument for the flare 
SOL19 80-06-21T02:00 (X 2.6) (upper limit) and then for the flare SOL1982-06-03T13:26 
(X8.0) jForrestetal.ll 19851) . The latter event provided the first convincing observations of 



Gamma Rays and Neutrons 



39 




06:40 06:44 06:48 06:52 06:56 

Universal Time (January 20, 2005) 



Fig. 3.1 Time profiles of the normalized flux density at 35 GHz (top) and normalized count rates of hard 
X-rays and y-rays (RHESSI, CORONAS-F/SONG) at different energies. The high ener gy y-rays observed b y 
SONG are pion-decay photons from primary protons at energies above 300 MeV (from Mass on et al.l2009T) . 

7T°-decay radiation in a solar flare. This event showed that the production of pions can oc- 
cur early in the impulsive phase defined by the production of hard X-rays near 100 keV. 
However, a significant portion of the high-energy emission was also observed well after the 
impulsive phase in an extended phase mostly observed at high energies. High-energy neu- 
trons were also detected for that event (see next section). Together with the development of 
models for ion acceleration and y-ray production, these observations led to the first deter- 
mination of high-energy ion distributions in flares jMurphvetaT]|l987h . Given the assumed 
spectral shape of energetic protons (either a Bessel function or a power-law spectrum), num- 
bers and spectra of energetic protons were estimated for both the impulsive and extended 
phases from the ratio of the differential pion decay radiation at 100 MeV to the 4.1-6.4 
MeV nuclear emission. It was found that the first phase was characterized by a significantly 
steeper proton spectrum than the extended phase. Pion-decay radiation was subsequently 
reported for other flares observed with SMM/GRS (Dunphy & Chupp 1991, 1992; see also 
Chupp & Dunphy 2000, for a review). 

After the end of the SMM mission in November 1989 and before the launch of CORONAS- 
F at the end of July 2001, there was no specifically solar-dedicated mission for high-energy 
radiation, but several satellites still provided y-ray data (GRANAT /PUEBUS, GAMMA- 1. 
and CGRO). 

Until now, around 20 events have been observed with significa nt pion production (see 
lLockwood et al.lll997l ; rMvagkova etafll2007l : IChupp & Rvar]|2009h . Some of the most re- 
cent events have been observed in a wide energy range by both RHESSI and CORONAS-F 
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Fig. 3.2 Left: 210 GHz submillimeter radio flux from KOSMA (Kolner Observatorium fur SubMillimeter 
Astronomie) and X-ray/y-ray fluxes from CORONAS-FISONG; right: imaging in radio at 210 GHz (red 
crosses), UV (background image from TRACE), X-rays and y-rays from RHESSI (blue and pink contours 
from Hurfo rd et al. 200 6) in phase B and C showing, respectively, electron and ion interaction sites (adapted 
from lTrottet et alj2008l) . 



(e.g jKuznetsov et al.ll200i : iTrottet et aHl2008l : iMasson et aljEoO^) (see, e.g., Figure [37Tb. 
For some of the events, high-energy emission has been observed for hours after the impul- 
sive phase of the flare, revealing that high energy ions mu s t be continuously accelerated 
on long timescal es in some flares (e.g jKanbach et al.lll993l : lRvan et alj[l994l ; lRvan|[2000l ; 
Ran k et all 1 200 lh . Quantitative analysis of a few of the events with significant pion pro- 
duction has been performed providing information on the ion energy spectrum at energies 
greater than 300 MeV and allowing a comparison of t his spectrum with the one d educed at 
lower energies from y-ray l i ne spectroscopy (see e.g., lAlexander et al.|[l994l ; iDunphv et al] 
ll999l : lKocharov et al.|[l994l 1 19981 : [Vilmer et al.ll2003l) . These comparisons generally show 
that the ion energy distribution does not have a simple power-law form from the y-ray line 
emitting energy domain (1-10 MeV) to the pion-producing energy domain (>300 MeV). 

In recent years, some of the events with significant pion-decay radi ation have been 
found to be associated with submillimeter (above 200 GHz) emissions ( Kaufmann et al. 
l2004l ; lLuthi et"aT]|2004l) . These observations have revealed the existence of a new emission 
component with intensities well above the extrapolation of the radio synchrotron spectrum 
seen at lower frequencies. Moreover, the radio spectrum increases with increasing frequency 
contrary to what is expected from an optically thin synchrotron spectrum. The origin of this 
new component is still under discussion. One possibility could be that this emission be pro- 
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duced by synch rotron emission from pion-dec ay positrons. This last process was described a 
long time ago l lLingenfelter & Ramatvlll967bl) and could be considered as a good candidate 
for the high frequency observations, given the association between events exhibiting pion- 
decay radiation and ev ents showing a spectral increase above 200 GHz iMyagkova et al.l 
120041 : ISilva et al.ll2007h . This is illustrated in Figure EOl from iTrottet et alJ ( I2008h which 
shows the observations of SOL2003-10-28T11:10 (X17.2) with KOSM^fat 210 GHz and 
with the SONG experiment on CORONAS-F at energies above 0.2 MeV. The onset of the 
impulsive component of the submillimeter component (phase B on the Figure) is clearly 
associated with the start of the radiation above 60 MeV, i.e., the start of >300 MeV/nucleon 
ion acceleration. The 210 GHz source size is compact (<10") in this phase and its location 
is co-spatial with the site of interacting ions (revealed by the 2.2 MeV line emission site) 
but not to the site of interacting electrons. In phase C, when no impulsive submillimeter 
emission and no strong emission above 60 MeV are observed, the submillimeter sources 
are quite different. The close correlation in time and space of the impulsive submillimeter 
emission and of the strong production of neutral pions thus suggests that synchrotron emis- 
sion from charged pion-decay positrons could be responsible for the submillimeter emis- 
sion. However the submillimeter flux predicted from the number of positrons derived from 
charged pion decay seems inadequate to account for observed fluxes. An alternative interpre- 
tation was given bv lSilva et all J2007h for SOL2003-1 1-02T09:55 (X3.9). It was concluded 
there that the most likely source of the submillimeter emission is gyrosynchrotron emission 
from non-thermal electrons but this requires high local magnetic fields, very large densities 
of accelerated electrons and very small sour ce sizes. Although many interpretations have 
been proposed in the literature (se e above and lKaufmann & R aulin 2006l lSakai & Nagasugil 
120071 : iFleishman & Kontarll2010h . the origin of this submillimeter component remains a 
challenging topic of discussion. 

Another challenging topic linked to the observations of pion decay radiation in a few 
flares is related to the origin of coronal hard X-ray so urces reported in, particularly, SOL2005- 
01-20T07:01 (X7.1). MacKinnon & Mallik (2010) have indeed shown that some coronal 
hard X-ray sources might be interpreted as photospheric, optical photons inverse Compton 
scattered to deka-keV energies by electrons, or indeed positrons in the 10-100 MeV en- 
ergy range. In particu lar, the coronal hard X-ray source reported in SOL2005-01-20T07:01 
(Kruckere t~al.ll2008bh could be explained in this way if just a few per cent of the secondary 
positrons implied by estimated fast ion numbers ( Masson et alj |2009) were present in the 
corona. Coronal hard X-ray sources in locations that seem too tenuous for a bremsstrahlung 
explanation (see iKrucker et al.l l2008al for a review) may thus offer a novel window on 
~GeV-energy ions. 



4 Neutron observations in space and on Earth 

Many reactions between energetic flare ions and ambient nuclei produce energetic neutrons 
as secondaries. These secondary neutrons undergo multiple, elastic scattering and either 
thermalize in the solar atmosphere or escape into interplanetary space. Measurements of the 
escaping neutrons augment knowledge of flare ion distributions gained from y-ray lines. In 
particular they fill a "diagnostic gap" between the de-excitation lines, dominated by ions in 
the 1-100 MeV energ y range, and the pion-decay radiation which needs ions of >300 MeV 
dLockwood et al.|[l997h . Neutrons in space give valuable information on the most energetic 
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flare ions, especially when measurements of pion decay y-rays are not available. They may 
also play a role in constraining the numbers of heavier fast ions, prolific producers of free 
neutrons via evaporation reactions. 

Having no electric charge, they arrive at the detector undeviated by active region mag- 
netic fields or interplanetary medium plasma turbulence. Free neutrons are unstable to f5- 
decay, with a rest-frame mean lifetime of 886 s, so many decay before they can be detected 
at Earth. 



4. 1 Neutron detections 

Energetic neutrons from flares have been detected both in space and on the ground. They 
give a signal in scintillators that may be distinguished by various means from photons or 
charged particles (pulse shape, comb ination of signals in multiple elements, etc.) (e.g.. 



Chu pp et al.lfl987l : iMurphv et alii 19991) . On the ground they have been detec ted by neutron 
monitors, and by specially designed neutron telescopes (e.g. JWatanabe et al]|2006b . 

The global neutron-monitor network has great value for studying extremes of solar par- 
ticle acceleration, offering detectors of a stopping power and effective area that would be 
impractical in space. They are located inside Earth' s magnetic field in a way that aids the dis- 
crimi nation of neutrons from other particle signals dUsoskin et alJl997l : lLockwood & Debrunnej 
Il999h . Neutron-monitor count-rate enhancements following solar fl ares are gene rally due to 
solar energetic particles interacting with Earth's atmosphere (e.g., Lopate 2006) but neu- 
trons genuinely of solar origin may be distinguished on grounds of timing and geographi- 
cal distribution. Neutron telescopes provide neutron spectral and angular res olution in the 
> 50 MeV energy range an d can discriminate neutrons from charged particles dMuraki et alJ 
ll992l ; lTsuchiva et al.l200ll) . They have been installed in seven high-altitude locations around 

the world. _ ^ 

A neutron of 150 MeV travels a distance of 1 AU in a neutron lifetime (MacKinnon 
120071) . Lesser energies are more likely to decay en route from the Sun while time dilatation 
extends the lifetime of neutrons at energies above 1 GeV. The neutron survival probability 
to Earth thus plummets below ~100 MeV. Nonetheless, solar flare neutrons have ar rived 
at Earth in d etectable numbers in the energy ranges 10-100 MeV dRvan et alH "l993) and 
50-360 MeV dMuraki et al.ll 19921) . Minimu m neutron energies f or detection on the ground 
are 50 MeV ( for neutron telescopes, see e.g., iMuraki et al.ll 19921) or 200-500 MeV (neutron 
monitors; see lDebrunner et al.ll 19891 : IShibatall 1994h . 

The first detection of flare neutrons was by t he SMM/GRS instrument, from the flare 
SOL1980-06-21T02:00 (X2.6) dChupp et al] 1982b . The flare SOL1982-06-03T13:26 (X8.0) 
produced energetic neutrons that were detected in space by SMM/GRS and on the ground 
using the Jungfraujoch and other neutron monitors (Chupp et all 1987b . Neutron-decay pro- 
tons were also detected in space from this flare IIEvenson et all 1983b . 

Gr ound -based detections of solar ne utrons have been recently reviewed bv lWatanabe et alJ 
d2005b and lValdes-Galicia et all d2009b . Since SOL1982-06-0313:26 (X8.0), ten more so- 
lar neutron event s have been clearly recognized in neutron monitor data: SOL 1990-05- 
24T21:45 (X9.3) dShea et al.ll99ll;ISmart et alll995l:lDebrunner et all 1991 1 1997b. SOL1991- 



03-22T22:45 (X9.4) (Pyle & S impson! 1 99 lb, SOL 1 99 1 -06-04T03 : 3 7 (X12.0) and SOL1991- 
06-06T03:37 (X 12.0) dMuraki et allll99a IStrurninskv et al.lll994b and several occasions 
durin g Cycle 23 dWatanabe et al.ll2003l : iBieber et all 20051 : IWatanabe et alll2006l : ISako et all 
2006). The second such event was observed by the IGY-type neutron monitor loc ated at Cli- 
max and several other stations in North America: SOL 1990-05-24T2 1:45 (X9.3) ( She a et all 
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ll99ll : ISmart et all 19951 : 15ebrunner et al.l 19911 19971) . In the SOL 1 99 1 -06-04T03 : 37 (X 1 2.0) 
event, solar neutron signals were recorded by three different detectors (neutron monitor, so- 
lar neutron telescope, and muon telescope) located at Mt. Norikura. 

In space, apart from the SMM/GRS instrument, energetic flare neutrons have been de- 
tected with the OSSE an d COMPTEL instruments on the CGRO mission dMurphv et alj 
1999 |; lDunphv et all 19991) as well as wi th the PHEBUS instrument aboard GRANAT dPebrunner 



1997 ) (see. e.g.. lchupp & Rvanll2009l for a revie w). Pulse-shape ana lysis in OSSE's Nal 
scintillators discriminated neutrons from photons (Mur phy et al.lll999l) . In COMPTEL, sig- 
nals from two arrays of detectors were combined to deduce arrival directions and energies 
for individual neutrons, via elastic scattering kinematics. Such imaging spectroscopy o f 10- 
100 MeV neutrons was possible for several flares in June 1991 (e.g. jRvan et al. 1993). 



4.2 Key features of neutron observations 

Solar neutron events are too few in number for st atistical studies. None theless detections 
so far share some common features summarized in lWatanabe et alj d2005l ). Detectable neu- 
trons have come from l arge flares, SOL2000-1 1-24T15: 13 (X2.3) being the smallest of these 
(Wata nabe et"aT]|2003r) . Of course, this could be observational bias rather than a physical 
feature of flares. No correlation is evident between the flare soft X-ray flux and the neutron 
emissivity. Detectable neutrons at Earth do not appear to come preferentially from limb or 
disk centre flares. 

Using, e.g., the y-ray de-excitation line flux time profile as a proxy for the rate of ion 
production at the Sun, neutron energy distributions have been deduced from neutron-monitor 
count-rate time profiles. Neutrons at the Sun are ge nerally found to have E~ s distributions, 
with 8 in the range 3-4 dWatanabe et~aT]|2003l . l200q) . Via model calculations using the code 
of iHua et al.l d2002h . the ion energy distribution at the Sun is then found to be roughly one 
power of E softer. Total emissivity at the Sun is in the range 4 x 10 28 -8 x 10 30 sr~' for 
neu trons in the range 0.05- 1.5 GeV. 

lLockwood et al.l ll 19971) and [Debrunner et al.l d 19971) emphasized that the de-excitation 
line flux does not always give the best proxy for the time development of neutron produc- 
tion, urging use instead of the pion-decay flux when this is measured. Neutron and pion 
production by energetic ions (>300 MeV) are closely related (see Section 1431 so this is 
not surprising. The flare SOL2005-09-07T09:52 (X17.0), for instance, was detected on the 
ground by both neutron telescopes and neutron monitors and found to produce neutrons over 
a more extended period than implied by the fl are impuls ive phase, as defined by hard X-ray 
emission. The class of extended y-ray flares (Rvan 2000) often show pion-decay y-radiation 
over tens of minutes. Similarly time-extended neutron emission would be a corollary, as was 
found u sing COMPTEL for SOL1991-06-09T04:24 (X10.0) and SOL1991-06-15T11:17 
(X12.0) ( lRvanetal.ll 19931) . 



4.3 Production of solar neutrons 

Ions are believed to be accelerated in the corona, losing their energy mostly in the chro- 
mosphere and photosphere and producing the bulk of their secondaries (photons, neutrons) 
there. The eventual escaping neutron velocity distribution is governed by the initial velocity 
distribution of neutrons, which in turn depends on the primary ion distribution throughout 
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the atmosphere, and the consequences of mutiple neutron elastic scattering. Thus a com- 
plete modeling of neutron production needs to include ion transport, neu tron production and 
scattering. Such a code was constructed bv lHua & Lingenfelted dl987bh and improved and 
updated by iHua et al. I J2002h . and t he consequences of varying loop and acceleration p a- 
rameters was described in detail by iMurphv et alj d2007h . The code of IHua et alJ J2002h is 
freely available and may be used to calculate 2.223 MeV line flux and neutron distributions 
expected at 1 AU under a wide variety of assumptions for both physical and acceleration 
parameters. 

Neutron-producing reacti o ns of fast ions ha ve been enumerated bv lHua & Lingenfeltej 
dl987bh : lMurphv et al.l ( fl987h : lHua et al]d2002h . The reaction p(p,n;z: + )p' has a threshold of 
270 MeV. It is a major contributor to total neutron yie lds for harder ion distributions (e.g., 
energy power law E~ s with s > 3: lMurphvet~aDl2007h . Thus neutron emission from flares 
will often be accompanied by pion-decay y-rays. 

In contrast are the cross-sections for collisions involving heavier nuclei. Neutrons in 
these cases can be produced via evaporation processes, without the creation of pions, so 
cross-sections often have much lower threshold energies, below 1 MeV/nucleon in some 
cases. Neutron measurements thus have some potential for constraining heavy ion acceler- 
ation in flares. Especially with softer assumed ion energy dis tributions, neutron y ields will 
depend sensitively on assumptions about fast ion abundances (Murphy et ~aDl2007n . 

Figure l4~Tl show s neutron distr ibutions produced by a population of fast ions, calculated 
using the code of H ua et al ] j2002h . A power-law ion energy distribution has been assumed, 
°=E~ S with s = 2, and fast ion abundances representative of "gradual" particle events in 
space. No magnetic field convergence or MHD pitch-angle scattering have been assumed 
in the corona so ions precipitate freely to the dense atmosphere. The two figures show the 
energy distribution of neutrons produced in the solar atmosphere, and the energy distribu- 
tion of neutrons that actually escape from the Sun. Above ~100 MeV the neutron energy 
distribution reflects the primary ion energy distrib ution; for assumed pow er laws the neutron 
distribution is about one power of energy harder dWatanabe et alJ [20051. At lower energies 
the neutron energy distribution is dominated by the relevant nuclear physics, particularly 
since many of the reactions are near threshold, and by neutron scattering in the atmosphere. 
Any factor that increases neutron production at greater depth (e.g., pitch-angle scattering 
in the strong limit which maxim izes precipitation - Section 12.2.6b decreases the escaping 
neutron flux ( IMurphv et al.ll2007l) . 



4.4 Neutron decay products 

Energetic protons fr om decaying flare ne utrons were detected in space fr om SOL1980- 
06-21 T02:00 (X2.6) fevenson et al]|l99Ch, SOL1982-06-03T1 3:26 (X8.0) dEvenson et all 
Il983h and SOL1984-04-25T01:40 (X13) dEvenson et al.lll990h . When the detector is not 
magnetically well connected to the flare site, neutron decay products may be clearly distin- 
guished from flare-associated fast particles. An energetic proton "precursor" in SOL1990- 
05-24T21:45 (X9.3), arriving at the GOES spacecraft before the bulk of the flare-associated 
fast i ons, has also been interpreted as a signature of decaying neutrons ( Kocharov et alj 
1 19961). Energetic ne utron-decay electrons were observed from SOL1980-06-21T02:00 (X2.6) 
dDroegeetal]|l996h . 

Neutron-decay protons carry off most of the energy of their pa rent particles and ar e 
straightforwardly detected, given an appropriately placed spacecraft dEvenson et al.| [T990). 
These advantages have to be weighed against the complication of accounting correctly for 
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Fig. 4.1 Angular-dependent energy spectra of the escaping neutrons at the Sun (solid curves) and of those 
surviving at 1 AU (dashed curves) at different cos 6 (zenith angle) compared to their production spectra at 
the Sun (dotted curves) in the case with no pitch-angle scattering of the ions (left) and of nearly s aturated 
pitch-angle scattering (right). See Section l2.2.6l for a discussion of scattering. (From Hua et al. 2002). 

proton propagation in the turbulent interplanetary medium dRuffololll99ll) . Clearly they 
could offer a useful window on flare ion acceleration but w e are aware of only one such 
claimed detection in Cycles 22 and 23 llKocharovetal.ll 19961) . and the argument in this case 
rested on fast proton arrival time rather than a lack of magnetic connection to the flare loca- 
tion. 



4.5 Low-energy neutrons 

Neutrons <100 MeV mostly decay long before they reach 1 AU. A neutron detector placed 
inside the orbit of Mercury would detect four to five orders of magnitude more neutrons in 
the 1-20 MeV energy range than a detector at 1 AU (in contrast to a detector of photons, 
which would gain roughly a single order of magnitude). There are se veral good scientifi c 
reasons to attempt neutron measurements in the inner heliosphere (e.g.. lVilmer et al.ll200lh . 
Much greater neutron fluxes will be measured from large flares by detectors on future inner- 
heliosphere space missions. As we mentioned above, heavy-ion populations may be usefully 
constrained, even though these do not produce narrow de-excitation lines. Constraints on 
ion acceleration in small flares and even the quiet Sun will be greatly strengthened. It should 
be emphasized that very recently, the MESSENGER neutron spectrometer made the first 
claimed detection of 1-8 MeV solar neutrons continuously produced for several hours, in 
good consistency with the previously observed e xtended production o f y-ray line emissions 
or pion-decay radiation (see previous sections) (Feldman e7aT]|2010h . However, local neu- 
tron production by solar energetic protons and a particles cannot be completely excluded in 
these observations. T hese consideration s motivate ongoing development of multiple scat- 
ter neutron detectors jPirard et al.ll2009h . whose angular resolution in particular is crucial 
to discriminate genuine solar neutrons from those produced locally by interactions of solar 
energetic particles in the spacecraft. 
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5 Imaging solar y-rays 
5.1 Introduction 

The spatial distribution of y-ray emission provides a channel of information on the acceler- 
ation and transport of energetic ions that is independent of spectroscopic analyses. As with 
electron bremsstrahlung, nuclear y-rays disclose not where particle acceleration occurs, but 
rather where accelerated particles interact with the ambient atmosphere. Nevertheless it pro- 
vides the only currently available method for placing accelerated ions in a solar magnetic 
context. 

Prior to the launch of RHESSI, the only information on the spatial distribution of so- 
lar flare y-ray emission originated from non- imaging observations of y -ray line emissions 
from flares believed to be behind the limb. IVestrand & Forrestl l ll993l) . for example, re- 
ported a strong neutron-capture line emission from a behind-the-limb flare. Such an ob- 
serva tion would imply that in addition to a compact source predicated in many models 
fe.g jHua & Lingenfeltejri987ah . the y-ray line emission also originated from a large scale 
source, probably diffuse, that was f ar from the flare i tself. Another event observed with 
strong prompt y-ray line emission bv iBaratetaD l l 19941) was also believed to be associated 
with a behind-the-limb flare. In this case, the absence of neutron capture line emission as 
well as the smooth time profile of the prompt y-ray line emission led to the conclusion that 
this emission was produced in the low corona at densities below 1 — 5 x 10 11 cm~ 3 . On the 
other hand, measurements of y-ray line emissions (time profiles with respect to HXR time 
profiles) have been interpreted as generally consistent with their production in a compact 
region located in the chromosphere at densities > 10 12 cm~ 3 (e.g., IChupplll996h . 

RHESSI has provided the first opportunity to locate the y-ray emission from solar flares 
directly. RHESSI y-ray imaging to date is based mainly on the narrow neutron-capture line 
at 2.223 MeV. This line is generated when fast neutrons, created by the interaction of ac- 
celerated protons and ions with the ambient atmosphere, are thermalized and subsequently 
captured by protons, producing deuterium and a narrow y-ray line at 2.223 MeV. Narrow 
because of the absence of Doppler broadening, the line is an optimal choice for imaging 
for four reasons: i) it is relatively intense; ii) since locally-generated fast neutrons cannot 
thermalize in the limited mass of the spacecraft, there is no locally-generated background 
line emission; iii) the line is narrow so that the continuum background is minimized; iv) 
the narrowness of the line also limits the contribution of the solar electron-bremsstrahlung 
continuum to the resulting images. The multip le elastic collisions r equired for neutron ther- 
malization require an average of ~100 s (e.g. jMurphv et al. 2007), which results in a char- 
acteristic time delay in this line compared to the prompt lines. Nevertheless, the correspond- 
ing average spatial offset between the primary interaction site and the capture site is only 
~500 km (R. Murphy, private communication). As a result, despite its two-step formation 
mechanism, a neutron-capture line image can faithfully indicate the site of the original nu- 
clear interaction to within an arcsecond. 

RHESSI imaging of nuclear y-ray emission is based on time- and energy-tagged counts 
from the thick rear-segments of RHESSI detectors 6 and 9 (35" and 183" resolution, respec- 
tively). These are the only detectors for which the grids are sufficiently thick (2 and 3 cm, 
respectively) to modulate the high energy photon flux. Rear segments of the detectors are 
shielded from the intense flux of low ene rgy X-rays, and so potential concerns about instru- 
mental live time or pulse pileup effects dSmith et alj ( 120021) can be discounted. On the other 
hand, since only two of the nine collimators (RMCs) are relevant to imaging and their thick 
grids block some of the flux, the y-ray throughput for such imaging is only about 14% of that 
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Table 5.1 RHESSI Events with Neutron Capture Line Imaging 



Parameter 


2002-07-23 


2003-10-28 


2003-10-29 


2003-11-02 


2005-01-20 


Time range 


00:27:20- 
00:43:20 


11:06:20- 
11:29:40 


20:43:00- 
21:00:00 


17:16:00- 
17:29:00 


06:44:00- 
06:56:00 


GOES class 


X4.8 


X17 


X10 


X8 


X7 


Energy range 
(keV) 


2223±5 


2223±5 


2223±5 


2223±5 


2223±8 


Total counts 
in RMC 6, 9 


336,240 


2280,1643 


463,313 


781,557 


961,518 


Relative visibility 
RMC 6 


1.61±0.45 


0.91±0.20 




0.76±0.26 


0.98±0.27 


Relative visibility 
RMC 9 


1.28±0.26 


0.82±0.10 


1.28±0.28 


0.84±0.15 


0.93±0.21 


Source size 
FWHM (arcsec) 
(2 o~ upper limit) 


<22 


<30 


<94 


<44 


<20 



for spectroscopy, and so the sensitivity for imaging is necessarily lower than that for spec- 
troscopy. The count data provide the st arting point for algori thms such as back-projection 
to reconstruct the image of the source dHurfordet all [2002). The count statistics and the 
relatively few spatial frequencies measured do not permit RHESSI to reproduce complex 
source morphologies. However, the data are well-suited to the quantitative characterization 
of simple sources and y-ray images for five flares have been obtained. 



5.2 Neutron-capture line imaging 

The complete set of images obtained to date is shown in Figure [BTTI with quantitative de- 
tails summarized in Table |5TI (Hurford et al. 2003, 2006a, 2006b). Rather than deal with 
each flare in turn, we summarize these observations from three perspectives: the number of 
distinct spatial components; their size; and their location. 

Considering first the number of distinct components, four of the five events presented a 
single, unresolved source whereas the fifth event, SOL2003-10-28T11:10 (X17.2), showed 
a double source that straddled an arcade of loops. This pattern does not necessarily imply, 
however, that double footpoint sources are an atypical feature of the y-ray emission. Two 
conditions must be met to image a double footpoint source: the separation of the footpoints 
must be sufficiently large that they can be resolved with the available 35 arcsec resolution; 
and there must be sufficient signal-to-noise that two sources (rather than one) can be confi- 
dently detected. Referring to the count statistics in Table 1 5.1 l and using the separation of the 
electron-bremsstrahlung components as a guide, SOL2003-10-28T11:10 (X17.2) was the 
only one of the five that satisfied both criteria and that event did display a double source. 
Had y-ray emission from the other events also been in the form of double sources, this 
could not have been seen. Therefore, based on the number of observed components, there is 
no evidence to suggest that the true neutron-capture line sources are predominantly single 
sources, a result that would be in contrast to typical hard X-ray double sources. 

A second perspective is provided by the estimates of the size of the unresolved compo- 
nents. Such estimates can be obtained from comparisons of the total observed line flux to 
that imaged in the unresolved source. The ratio of imaged flux to total flux, called the rela- 
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Fig. 5.1 Respective locations of y-ray sources at 2.223 MeV and of hard X-ray sources observed with RHESSI 
for different flares. The SOL2002-07-23T00:35 (X4.8) event; the circles represent the 1 c errors for the 300- 
500 keV (black), 700-1400 keV (red) and 2218-2228 keV (pink) maps made with identical parameters. The 
white contours show the high resolution 50-100 keV map with 3 arcsec res olution. The background is a 
TRACE image showing the post-flare loops (adapted from lHurford et alj|2003l) . The SOL2003-10-28T11:10 
(X17.2) event; overlay of contours at 50%, 70%, and 90% of 35"y-ray images made with detectors 6 and 9 
on a TRACE image. The red plus signs indicate 200-300 keV footpoint locations at different times. The 
SOL2003-10-29T20:49 (X10.0) event; fine red contours: (50%-90%) of the 200-300 keV high resolution 
maps superposed on a TRACE image. The purple line indicates the motions of the 200-300 keV footpoints. 
Thick red and blue circles: centroid locations with 1 a error for detector 9 imaging at 200-300 and 2218- 
2228 keV, respectively. The SOL2003-1 1-02T17:25 (X8.3) event; overlay of 200-300 keV (red contours) 
and 2218-2228 keV (blue contours) sources on a postflare SOHO/EYT image. The dashed and thin solid 
red contours (50%) show the result of 200-300 keV mapping with 7" and 35" respectively. The thin blue 
contour (50%) is the corresponding map at 2218-2228 keV made wi th detectors 6 and 9. T he thicker lo" error 
circles show the detector 6 source centroid location (adapted from Hurford et al. 2006b) for the SOL2005- 
01-20T07:01 (X7.1) event; contours of the 250-500 keV emission overlaid on a TRACE image and centroids 
(in red and blue) of the map performed in the 250-500 keV and 2215-2231 keV range. 
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tive visibility, can be used to put upper limits on the size of an unresolved source. Ignoring 
the modulation, the total number of observed counts in the line provides a direct indication 
of the total, spatially-integrated flux in the narrow energy range associated with the line. 
Such counts are due to a combination of real solar neutron-capture line events, continuum 
solar emission and background. Since the number of such counts that are due to solar con- 
tinuum or background continuum can be reliably estimated from energy ranges just above 
2.223 MeV, an accurate determination of the spatially integrated solar neutron-capture line 
flux can be made, independent of background models. The imaged flux does not contain 
any unmodulated background, but could include a contribution from the solar continuum. 
Although this contribution could be estimated by imaging in a narrow energy range above 
2.223 MeV, it is easy to demonstrate that, thanks to the narrow energy range required for 
imaging the neutron-capture line (10-16 keV), its contribution is negligible. Note that since 
the relative visibilities are estimated by comparing the imaged to total counts observed with 
a given rotating modulation collimator, the measurement of relative visibilities is indepen- 
dent of potential errors in estimating the detector efficiency. 

The observed relative visibilities (Table 15.1b show that in all cases the flux in the un- 
resolved neutron-capture line source(s) is consistent with the spatially-integrated flux. This 
implies that there is no evidence for a large-scale diffuse component. By assuming a Gaus- 
sian source model, upper limits to the FWHM source diameters are obtained that correspond 
to a few tens of arcseconds. 

A third perspective on the neutron-capture line sources is provided by their location. An 
unanticipated result was that in some, if not all, cases, the neutron-capture line sources were 
not co-spatial with the corresponding electron-bremsstrahlung sources. Since the character 
and statistical significance of the displacements differed from event to event, we first con- 
sider them individually. For SOL2003-10-28T1 1 : 10 (X17.2), double footpoint y-ray sources 
straddled the same loop arcade as the electron sources. However, they were displaced along 
the loop arcade by 14" and 17" ±5" respectively (Figure [5TTT l. For the other events, the loca- 
tion of the single unresolved sources are expressed as a centroid location with appropriate 
statistical errors. For SOL2003-07-23T00:35 (X4.8), the neutron-capture line source was 
displaced by 25"±5" from the corresponding electron source. Although it was not located 
near any obvious EUV or Ha line emission, its location was consistent with the footpoint 
of a post-flare loop. For SOL2005-01-20T07:01 (X7.1) the neutron source was clearly asso- 
ciated with one of the two electron-bremsstrahlung footpoints. For SOL2003-11-02T17:25 
(X8.3) the location also appeared to be preferentially associated with one of the two electron- 
bremsstrahlung footpoints, although in this case the displacement of the 2.2 MeV emission 
compared to the centroid of the electron-bremsstrahlung source was significant only at the 
2a level. For SOL2003-10-29T20:49 (X10.0), the measured separation of the correspond- 
ing centroids was 13"±17", so a separation, if any, could not be established. In summary, 
statistically significant displacements were observed in three of the five events, marginally 
in a fourth but not seen in the fifth and statistically least significant event. 

5.3 Other y-ray imaging 

In addition to the neutron-capture line, imaging of the y-ray continuu m in the 200-800 keV 
range dominated by electro n bremsstrahlung has also been reported (Hur fordetai1l2005L 
2006a; Krucker et al. 2008b). At higher energies, where nuclear emission is more important, 
only RHESSI detector 9 (183") was able to modulate effectively with the events studied. 
The tentative explanation for the reduced efficiency of detector 6 modulation for imaging 
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in a broad energy band at high energies is that photons passing through front grid slits into 
rear grid slats can scatter into the detector and be counted as part of the modulated flux. 
The modulation of such events, however, is 180° out of phase with that of the valid (front 
slit-to-rear slit) events and so tends to reduce the effective modulation efficiency. This issue 
is a factor for detector 6 (but not detector 9) since the detector 6 grids are not thick enough 
to suppress the scattered photons. This issue is not relevant for neutron-capture line imaging 
since photons that are scattered by the front grid scatter away from the narrow window used 
for spectral analysis. 

Nevertheless, the centroid locations obtained for the SOL2005-01-20T07:01 (X7.1) y- 
ray continuum with detector 9 in four bands between 600 keV and 8 MeV were of interest 
particularly because of their energy dependence. At the lowest energies in this band, their 
centroid locations matched that of the centroid of the X-ray electron-bremsstrahlung source. 
However, at higher energies where the nuclear contribution became more important, the 
y-ray c ontinuum location trended toward the neutron-capture line centroid l lHurford et al.l 
l2006ah . Such a pattern is consistent with the displacement between the sites of electron- 
bremsstrahlung collisions and nuclear interactions. 

Imaging in a narrow energy range surrounding the 511 keV line was also attempted 
for SOL2005-01-20. Images obtained independently from the front and rear segments of 
detector 9 were found to be self-consistent in terms of strength and intensity. However, 
modeling suggests that although a narrow 500-520 keV energy range was used, the imaged 
fluence was dominated by the electron-bremsstrahlung continuum. 

5.4 Discussion 

One well-established result of the neutron-capture line mapping was that the flux for all five 
events was consistent with coming from compact source component(s) in the flaring active 
region. This would suggest that the flare-associated y-ray emission was associated with flare- 
accelerated ions, and not with ions that are shock- accelerated at high altitudes. If that were 
the case, the ions might be expected to interact with ambient material in a more diffuse 
pattern, perhaps at some distance from the flare. This shows that the y-ray emission does not 
result from the interaction of solar energetic particles with the solar wind as alternatively 
proposed bv lKahler & Ragotl d2008h . 

Given the resolution and sensitivity of the imaging and the limited sample of events, 
there is no evidence that the morphologies of nuclear y-ray sources (e.g., double vs. single 
footpoints) differ systematically from that of hard X-ray sources. 

The most striking result, of course, is the differences in the locations of the electron- 
bremsstrahlung and neutron-capture line sources. Hypotheses that attempt to account for 
such differences generally fall into three catagories: i) the displacement is an instrumental 
artifact; ii) the displacement has a solar, but relatively mundane origin; or iii) the displace- 
ment reflects differences in acceleration and/or transport processes for flare-accelerated ions 
and electrons. 

We first consider the possibility that the displacements are an artifact of the imaging 
process. There are three independent arguments indicating that this is highly unlikely. First, 
the measurements of the locations of the sources at two different energy regimes were based 
on analyses that differed only in the energy of the counts upon which the imaging is based. 
The same time interval, detector grid optics, aspect solution, imaging algorithm and imaging 
parameters were used for both so that the potential for systematic co-location errors can 
be discounted. Second, if the relative displacements had a systematic instrumental origin 
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they would likely be with respect to the only axis of symmetry associated with the imaging 
process, namely the line joining the direction of the rotation axis to that of the flare. However, 
the displacements were neither systematically parallel to or perpendicular to this direction. 
Third, in all but SOL2002-07-23T00:35, the direction of the displacement, if any, was in a 
magnetically significant solar direction, as opposed to some randomly-chosen direction. As 
result of these three considerations, it is considered highly unlikely that the displacements 
have an instrumental origin. 

A second class of potential explanations for the displacements is that they have a solar, 
but relatively mundane explanation. One such explanation might be based on the average 
distance between the locations of the original nuclear interaction and that of the neutron- 
capture process. However, as mentioned, such a displacement is expected to be less than 
Another possibility is that the displacement is related to the time-delay associated with 
thermalization. This scenario notes that the neutron-capture line image reflects the location 
of nuclear interactions that occurred ~100 seconds earlier. To investigate this hypothesis, for 
each event the neutron-capture line images were compared both to cotemporaneous electron- 
bremsstrahlung images and to electron-bremsstrahlung images made with integration times 
shifted by 100 s. No significant differences were observed. 

These considerations therefore support a solar-physics-based interpretation of the dis- 
placement, based on a difference in acceleration and/or propagation between accelerated 
ions and electrons. Several such interpretations have been put forward. lEmslieetalJte004bh 
discussed a stochastic acceleration mechanism whereby ions a re preferential ly accelerated 
in larger loops than electrons. Imaging for SOL2002-07-23 dHurford etal.ll2003l) would 
support this interpretation. The interpretation of time delays between y-ray lines and hard 
X-rays by trap-plus-precipitation models for the same flare ( IDauphin & Vil mer 2007) also 
confirm the fact that ions are injected in longer loops than electrons. 

However, such an explanation would not appear to be appropriate for the SOL2003- 
10-28T1 1:10 (X17.2) where the footpoint separations for t he electrons and ions w ere sim- 
ilar. Gradient and curvature drifts were also considered by iHurford et al.l j2006bl) . but the 
predicted displacements between th e ions and electrons were tw o orders of magnitude too 
small. Another possibility, raised bv lLitvinenko & Somovl dl993h . was that electric-field ac- 
celera tion would lead to separation of particles with opposite charge dZharkova & Gordovskwl 
2004). Its potential applicability in these cases has not yet been pursued. 

In two of the five events (SOL2005-01-20T17:25 and SOL2003-11-02T17:25) the ions 
were preferentially associated with one of the two electron bremsstrahlung footpoints. This 
raises the possibility that at least for these two events, the different footpoint weighting 
between the two species might be attributed to differences between electron and ion prop- 
agation in asymmetrical fields. This would be a result of the significant differences in the 
scattering properties and ionization-loss stopping distances between the accelerated elec- 
trons and ions. Modeling to support such a hypothesis has not been carried out, however. 

With only five imaged events, it is too early to establish the systematics of the dis- 
placements. It is also apparent that RHESSPs 35" resolution is insufficient to resolve the 
individual y-ray sources. Therefore, it is to be hoped that RHESSI will observe additional 
events in the coming maximum so that their systematic properties can be further explored 
and that new instrumentation will become available with improved angular resolution and 
sensitivity in the y-ray regime. 
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6 Conclusions and future perspectives 

Two particular observational advances characterize the RHESSI era: spectral resolution ad- 
equate to reveal line shapes and to better constrain the line fluences; and y-ray imaging. 
The former provides in principle previously inaccessible information on proton and a- 
particle energy and angular distributions. The latter has revealed the sites of ion interac- 
tion and neutron production and most dramatically shown that they do not coincide with 
the HXR-identified locations of electron interaction. Also, it has not confirmed the SMM- 
era suggestion that 2.223 MeV line emission might be produced over an extended region 
l lVestrand & Forrestll 1993b . The observations are also not consistent with a production of 
y-rays through the inte raction of solar energetic protons with the solar wind as examined by 
lKahler&Ragol(l200 8). We may look forward to further advances, both observational (e.g., 
with GRIPS uShih et alj |2009a). and theoretical, over the next few years. 

At the same time as RHESSI as well as INTEGRAL/SPI observations have brought pi- 
oneering observations of y-ray spectra at high resolution, it was also realized that the de- 
termination of all the parameters which characterize the y-ray spectrum (target abundances 
and ionization states, energy and angular distributions of energetic ions, relative abundances 
of accelerated a particles and heavier ions, transport and energy losses of energetic ions be- 
tween the acceleration region and interaction region if not co-spatial, etc.) is still very chal- 
lenging. In particular, the different parameters cannot be determined independently from the 
assumptions made on others: both ambient abundances and parameters of energetic particles 
must be deduced together even in the "simple" analysis of a y-ray spectrum. Combining the 
information from narrow line shapes and line fluence ratios can, however, help to disentan- 
gle the effects of the different parameters. The use of lines purely produced by a particles 
(in particular the line at 0.339 MeV formed by the interaction of a particles with iron) would 
also considerably improve the analysis of y-ray spectra by providing cleaner constraints on 
the distribution of a particles. A few attempts had been made at the time of SMM, but no 
further studies have been published since then based on spectra obtained at higher resolu- 
tion. This is also the case for the lines produced by energetic 3 He on ambient 16 around 
1 MeV. Line shapes and widths of the 511 keV annihilation line and positronium continuum 
have been measured at high spectral resolution leading for a few flares to unexpected com- 
bination of chromospheric densities and transition region temperatures for the annihilation 
region. Taking this into account in the future, the production and annihilation of positrons 
over a broad range of heights and thus of conditions in the medium is probably the next step 
to undertake to better understand the observations of shapes and widths of the annihilation 
line and continuum. 

Heavier (Z > 2) accelerated ion species are harder to diagnose but they may non ethe- 
less be significant in overall flare energetics dMurphv et all 19971 : lEmslie et alj|2004 sb . The 
detailed treatment now availa ble of the unresolved component of nuclear y-rays should bear 
significantly on this question jMurphv et al.l2 009). Heavier species are prolific producers of 
free neutrons via evaporation-type reactions so instruments planned to measure, on future 
inner -heliosphere missi ons, neutrons in the 1-10 MeV range not yet explored systematically 
will yield new constraints. 

Do all flares accelerate ions? It seems unlikely that they do not, but all flares do not 
produce detectable y-ray line s. RHESSI' s eight years of operation have accumulated a data 
set bearing on this question (Shihetal. 2009b). A good correlation is found between the 
energy content of electrons above 300 keV and protons above 30 MeV suggesting a strong 
correlation between the acceleration of high-energy electrons and y-ray line producing ions. 
A larger variation is found between the energy content of electrons above 20 keV and the 
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ion energy content. Some episodes in flares or some specific flares exhibit indeed higher 
relative numbers of electrons with respect to ions. Based on the RHESSI data set, it was also 
found that there is a strong link between the production of energetic ions and sub-relativistic 
electrons, but only when a threshold in the production of ions is attained. In conclusion, 
while the acceleration of protons above 30 MeV is related to the acceleration of relativistic 
electrons, the acceleration of sub-relativistic electrons is linked to the acceleration of these 
energetic protons and relativistic electrons, but only when a threshold of energetic ions and 
electrons is reached. These conclusions must be kept in mind when addressing acceleration 
mechanisms in flares. The content in low energy ions from flares is mostly unknown and no 
strong constraints have been derived on the production of low energy protons in solar flares. 
Attempts to constrain fluxes from low energy proton capture lines could be a way to address 
this issue. Another promising way could be to detect suprathermal ions through charge ex- 
change with the ambien t hydrogen (see the first detection of energetic neutral hydrogen 
atoms in a solar flare by Mewaldt et al. 2 009h . 

Apart from CORONAS-F/SONG (e.g jKuznetsov et alj|2006h . solar observations at the 
highest y-ray energies have been lacking in recent years. This is one of the only ways (with 
the detection of neutrons) to address the issue of the highest energetic particles accelerated 
in solar flares and to determine the shape of the energetic ion spectrum over several decades 
of energy range. The global network of neutron detectors thus plays at the present time 
a particularly important role, yielding other wise unobtainable i nformation on ion energy 
distributions in the GeV energy range (e.g. jMasson et al ]|200i). As the solar cycle picks 
up, Fermi observations of the high-energy continuum should become critical to probing the 
highest ion energies attained in flare s. Further observa tions of flares at submillimeter (and 
even shorter) wavelengths (see, e.g. jKlein et alJ 20 ^) as well as the s earch for X-ray in- 



verse Compton radiation in coronal sources ^MacKinnon & Mallikll2010h could additionally 
provide information on electrons and positrons produced in the 10-100 MeV energy range 
in flares. 

Recent years have seen a re-examination o f the thick-target interpretation of flare HXRs 
dFletcher & Hudson! [20081 : brown et al.ll2009h . Since HXR-emitting electrons are believed 
to embody a large fraction of the total energy released in the flare, this is a question that goes 
beyond interpretation of radiation signatures to the heart of the flare energy release process. 
The existing thick-target models assume that electrons radiate negligibly in the (coronal) ac- 
celeration region, and that they slow down collisionally in the (thick-target, chromospheric) 
interaction region. If acceleration and interaction regions were one and the same, fewer 
electrons would be needed to account for HXRs. A similar reconsideration of y-ray source 
regions may prove necessary. 
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